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INTRODUCTION 


Sleep  has  been  studied  from  a  very  large  number  of  points 

of  view;  from  studies  involving  sleep  deprivation,  examination 

of  sleep  hemodynamics,  studies  examining  penile  erection  during 

(51 ) 

sleep,  schizophrenia,  gastric  motility  and  circadian  rhythms 

and  even  biochemical  and  seasonal  influences  on  sleep.  The 

vastness  of  research  conducted  on  sleep  can  probably  best  be 

perceived  after  reviewing  some  of  the  4,300  references  in  the 

(35  ^ 

bibliography  prepared  by  Nathaniel  Kleitman  '  entitled, 

"Sleep  and  Wakefulness". 

}  Our  own  investigation  into  the  relationship  between  sleep 
depth  and  instantaneous  beat-by-beat  heart  rate  has  been  carried 
out  intermittently  over  the  past  ten  years^  '  '  '  .  Many 

properties  of  the  relationship  between  heart  beat  patterns  and 
the  quantity  and  quality  of  sleep  have  been  investigated  in  our 
laboratory.  \ 

The  technical  advances  in  small  computers  has  created 

/ 

an  upsurge  ^n  research  to  find  an  automated  method  of  deter¬ 
mining  sleep  quantity  and  quality.  While  the  hardware 
capabilities  of  microprocessors  have  reached  the  stage  of 
feasibility  for  automated  sleep  staging,  software  development 
has  not  kept  pace.  Thus,  an  automated  on-line  sleep  evaluation 
device  still  remains  a  desirable  yet  unattained  reality. 

^-Eor  the  past  decade,  a  number  of  projects  have  been  con¬ 
ducted  in  the  Biomedical  Engineering  Program  here  at  the 
University  of  Texas  regarding  the  development  of  automated 


;  sleep  scoring  software.  Initially  our  studies  were  prompted 
by  the  advent  of  the  space  age  where  it  became  important  to 
devise  a  method  of  determining  levels  of  alertness  of  astro- 

Nj>/ 

nauts  during  prolonged  space  flight.  Specifically,  these 
projects  involved  bandwidth  reduction  of  sleep,  extraction 
of  sleep  information  from  heart  rate  data,  analysis  of  sleep 
cycles,  detection  of  REM,  1  sleep  stage  of  eye  movement  from 
beat-by-beat  heart  rate,  classification  of  sleep  into  awake, 
REM- 1 ,  and  Stage  2,  3  and  4  using  the  "beatquency"  domain. 

All  of  these  projects  were  undertaken  with  a  common  goal  in 
mind.  That  is,  the  development  of  an  automated  process  by 
which  rapid,  inexpensive  determinations  of  levels  of  alertness 
could  be  determined  accurately  using  an  easily  obtained 
physiologic  parameter.  The  physiologic  parameter  chosen  in  our 
studies  was  the  beat-by-beat  heart  rate. A  Our  investigations 
have  produced  some  insight  into  this  probjlem,  but  we  have  not 


BACKGROUND 


(17  26) 

We  do  not  really  know  what  sleep  is'  '  ;  however,  we 

all  sleep  approximately  one-third  of  our  lives  away.  We  know 
that  if  we  do  not  sleep,  we  function  poorly,  and  if  the  quantity 
and  quality  of  our  sleep  is  signif icantlv  altered  we  do  not 
function  with  our  usual  capabilities ^27 ' 31 ' 55^ . 

With  the  advent  of  the  polygraph  and  sleep  staging,  a 
method  of  quantitation  and  gradation  of  sleep  became  possible. 
While  it  would  seem  true  that  some  insight  into  the  nature  and 
function  of  sleep  might  result  from  this  quantification,  sleep 
staging  has  merely  raised  additional  unanswered  questions  about 
sleep  quality. 

Electroencephalographic  (EEG)  patterns  obtained  from  poly¬ 
graph  recordings  exhibit  definitive  characteristics  in  both 
awake  and  sleeping  stages.  Descriptions  of  these  patterns  employ 
objective  terms  such  as  alpha,  beta,  delta  and  theta  waves,  low 
voltage  fast  electroencephalographic  waves,  high  voltage  slow 
waves,  spindle  activity,  K-complexes  and  the  like;  however,  EEG, 
like  sleep  itself,  exists  as  a  continuous  process  rather  than 
as  a  sequence  of  four  discrete  states. 

With  the  discovery  of  ocular  motility  during  sleep ^ ' 5 ' 19 ' ^ ^ 
another  measure  became  available  for  objective  observation.  The 
electro-oculogram  (EOG)  measures  electrical  activity  caused  by 
eye  movements.  The  correlation  of  concurrent  stage  1  EEG  sleep 
and  EOG  patterns  with  dreaming  led  to  the  definition  of  REM 
(rapid  eye  movement)  sleep. 


The  classical  technique  for  determining  sleep  levels  or 

sleep  stage  involves  the  clinical  interpretation  of  concomitant 

changes  in  the  EEG  and  EOG  patterns  of  polygraph  traces  on  a 

continuous  all  night  strip  chart  record.  With  these  findings, 

researchers  have  developed  standard  procedures  for  "scoring" 

EEG  patterns  into  different  depths  of  sleep.  This  scoring 

process  is  a  standard,  somewhat  tedious  and  time  consuming 

process.  Nevertheless,  this  method  has  provided  us  with 

quantitative  techniques  for  research  and  has  allowed  us  to 

perform  extensive  analyses  of  sleep  and  its  related  phenomena. 

While  each  sleep  researcher  tends  to  develop  his  own  method 

of  quantitating  sleep  depth,  sleep  is  typically  divided  into 

four  primary  stages:  Stages  awake  (or  zero) ,  I,  II,  III  and  IV, 

with  Stage  IV  representing  the  deepest  sleep .  Stage  REM 

sleep  (labeled  Stage  V)  produces  the  same  EEG  patterns  as  Stage 

I,  but  during  REM  sleep  there  is  intermittent  rapid  conjugate 

movement  of  the  eyes.  Typical  patterns  for  EEG  and  EOG  seen 

( 35 ) 

during  these  stages  are  shown  in  Figure  1  .A  graphical 

representation  of  staging  throughout  the  entire  night  of  sleep 

( 39 ) 

is  shown  in  Figure  2  .  This  figure  depicts  the  sleep  subject 

progressing  from  Stage  awake  or  Stage  IV  and  then  returning 
rapidly  to  Stage  I,  REM,  which  is  Stages  I  and  REM  combined. 
This  pattern  is  typical  of  a  "normal"  night  of  sleep.  The 
horizontal  bars  above  Stage  I,  REM,  indicate  the  occurrence  of 
rapid  conjugate  eye  movements,  or  Stage  REM  (V).  It  can  also 
be  seen  that  the  subject  makes  a  similar  progression  from  light 


w*  ^  +  S’ 


sleep  to  deep  sleep  throughout  the  night.  This  "cycling"  occurs 
approximately  every  60-100  minutes  and  has  been  described  in 
numerous  reports  .We  have  defined  the 

cycles  as  beginning  at  the  onset  of  the  REM  period  (horizontal 
bars)  and  ending  at  the  onset  of  the  next  REM  period,  except 
for  Cycle  0,  which  typically  begins  with  Stage  I  where  there 
are  no  rapid  conjugate  eye  movements.  The  sleep  subjects' 
record  in  produced  five  complete  sleep  cycles  during  this  night 
of  sleep.  Usually,  with  each  sleep  cycle  the  "trough"  of  the 
cycle  gets  more  shallow  as  the  night  progresses,  which  is  also 
typical  of  the  normal  night  of  sleep^^'^^.  In  general. 


1 


(16) 


progression  from  wakefulness  to  Stage  IV  at  the  beginning 
Cycle  0  is  a  continuum  of  change,  while  reversion  to  the 
lighter  stages  at  the  end  of  the  cycle  are  quite  abrupt' 

Stages  III  and  IV  are  only  occasionally  reached  during  sleep 
cycles  following  the  second.  Another  characteristic  of  sleep 
cycling  is  that  the  length  of  the  REM  periods  increases  as 
the  night  of  sleep  progresses,  while  non-REM  durations 
decrease'4'14'25’. 

A  considerable  amount  of  research  is  currently  in  progress 
dedicated  to  a  better  understanding  of  the  relationship  of 
sleep  patterns  and  heart  rate  patterns  in  infants  who  are  apt 
to  experience  unexplained  crib  death ' ^ ^ ^ t 47 , 56 ) .  jn 
addition,  it  is  well  known  that  sleep  patterns  of  man  differ 
from  those  of  other  mammals '  ^  ) ;  however,  sleep  in  man, 


while  variable  from  night  to  night,  is  relatively  consistent  in 


terms  of  the  pattern  of  sleep  cycling  within  one  subject  and 

from  one  subject  to  another.  Age  plays  an  important  role  in  the 

(21) 

characteristics  of  sleep  patterns  . 

It  should  be  remembered  that  although  sleep  is  described  in 
quantitative  measures  or  "stages,"  sleep  patterns  are  associated 
with  a  continual  process.  This  means  that  sleep  levels  do  not 
always  lend  themselves  to  discrete  classifications.  In  other 
words,  there  are  fluctuations  of  depth  within  each  stage  as  well 
as  between  stages ^ 16 ' 63 ' . 

The  measurement  of  a  variety  of  physiologic  parameters 
associated  with  autonomic  activity  are  known  to  be  associated 


with  changes  in  sleep  patterns  15, 17, 27, 32, 35, 43, 52, 53, 54, 63) 


In  fact,  the  depth  of  sleep  is  said  to  take  on  different  meanings 
according  to  the  variables  used  to  define  it.  The  determination 
of  depth  should  take  into  account  the  simultaneous  fluctuations 
of  heart  rate,  blood  pressure,  respiration  and  even  skin 
potential.  The  measurement  of  gastric  motility  during  sleep, 
metabolic  rate,  penile  erections,  peripheral  muscle  activity 
and  pupillary  reactions  have  also  been  regarded  as  autonomic 


correlates  of  sleep.  Before  the  advent  of  the  polygraph,  auto- 

(63) 


nomic  variables  were  the  only  means  of  describing  sleep 
Summarizing  as  best  is  possible  the  commonality  in  these  auto¬ 
nomic  correlates  of  sleep  depth  research  are: 

1)  In  general,  the  trend  is  for  heart  rate,  blood 


pressure  and  respiratory  rate  to  decrease  with  increasing 
depth  of  sleep. 


2)  In  general,  heart  rate,  blood  pressure  and 


respiration  increase  their  variability  during  REM  sleep  and 
decrease  their  variability  with  increasing  depth  of  sleep. 

3)  Gross  body  movements  exhibit  cyclic  variations 
closely  associated  with  sleep  cycles. 


4)  Skin  potentials  generally  increase  with  the  onset 
of  deep  sleep  and  are  markedly  reduced  during  REM  si*-  'p. 

Concerning  heart  rate,  in  general  it  is  found  at  the 
heart  slows  from  two  to  twelve  beats  per  minute  wit  the  onset 


of  sleep 


(52) 


And  on  the  average,  the  range  of  hea.  *  .ate  is 


lower  during  sleep  than  during  waking.  Virtually  all  researchers 
(1,3,10,11,52,54,63) 


agree 


that  the  variability  of  heart  rate  is  the  most 


pronounced  feature  when  considering  different  depths  of  sleep, 
since  average  heart  rates  overlap  for  most  stages  of  sleep. 
Fluctuations  in  the  heart  rate  are  more  erratic  during  light 
sleep  as  in  Stages  I,  II,  and  REM  and  the  fluctuations  are 
slow  and  rhythmic  in  deeper  sleep  stages  such  as  Stages  III 
and  IV.  REM  sleep  is  associated  with  the  most  outstanding  and 
violent  variations  in  heart  rates;  yet,  paradoxically,  it  also 
contains  periods  of  almost  complete  heart  rate  quiescence (1° ' 52) . 

Both  blood  pressure  and  respiratory  rate(11)  follow  a 
pattern  and  trend  similar  to  heart  rate  during  the  sleep 
cycling  process.  The  concensus  view  is  that  blood  pressure  and 
respiration  show  high  degrees  of  variability  during  REM  sleep 
and  become  more  stable  in  deeper  stages  of  sleep  ^ "  »15'25'33,52,63) 

On  the  average,  blood  pressure  and  respiratory  rate  during  REM 
sleep  are  significantly  higher.  Blood  pressure  tends  to  reach 
its  minimum  during  the  first  third  of  a  night  of  sleep  and  is 


probably  associated  with  the  preponderance  of  Stages  III  and  IV, 
which  occur  almost  exclusively  during  this  period  of  the  night. 
Respiratory  changes  such  as  increased  respiratory  rate  and  de¬ 
creased  amplitude  have  been  correlated  with  REM  activity  (periods 

when  the  eyes  are  actually  moving)  v  '  '  .  Body  movements  occur 
periodically  throughout  the  night  particularly  associated  with 
the  deeper  sleep  at  the  end  of  each  sleep  cycle.  Typically, 
body  movements  are  seen  to  occur  during  Stage  IV  sleep  just 
prior  to  REM^^'^^.  The  REM  period  is  characterized  by  no 
peripheral  activity,  but  movement  begins  again  immediately 
after  the  REM  and  then  returns  toward  quiescence  with  deepenincr 
sleep.  The  almost  complete  body  paralysis  during  REM  periods 
exists  in  sharp  contrast  to  the  violent  autonomic  activity 
associated  with  REM  periods  and  is  probably  necessary  to 
avoid  destructive  behavior  and/or  self-injury  during  REM  sleep. 

Changes  in  functions  mediated  by  autonomic  information 
such  as  heart  rate,  bLood  pressure  and  respiratory  rate  are 
thought  to  follow  patterns  mediated  through  parasympathetic 
and  sympathetic  neural  control  In  some  cases,  these 

variables  appear  to  be  interdependent  upon  one  another.  For 
example.  Bond  et  al.^^,  found  that  heart  rate  variability 
during  Stages  I  and  II  was  rhythmic  and  roughly  associated 
with  respiration.  During  deeper  sleep,  as  in  Stages  III 

and  IV,  there  appeared  to  be  an  even  intercorrelation  and 

synchrony  between  heart  rate  patterns  and  respiratory  activity. 

However,  during  REM  periods  there  was  total  disassociation 

( 52 ) 

between  heart  rate  and  respiration.  Snyder,  et  al.  ,  also 

reported  that  heart  rates,  blood  pressure  and  respiration 

"oscillated  quite  regularly  around  the  relatively  stable  base- 

11 


line"  during  Stages  II  and  III/IV  combined,  while  Stages  I  and 
REM  combined  featured  wide  and  erratic  fluctuations.  Snyder 
reported  many  instances  where  changes  in  all  three  variables 
occurred  together  and  in  other  instances  completely  independent 
change.  They  suggested  that  simultaneous  changes  of  these 
variables  are  probably  secondary  to  changes  in  respiration. 

Buast ,  et  al.  performed  detailed  experiments  on  parasympathetic 
and  sympathetic  influences  on  heart  rate  during  sleep  using  the 
cat.  He  concluded  that  phasic  changes  in  the  heart  rate  during 
sleep  and  wakefulness  are  decisively  influenced  by  autonomic 
events.  in  general,  he  found  that  the  fall  in  heart  rate  from 
wakefulness  to  synchronized  (Stages  II,  III,  iv  in  the  human) 
sleep  is  due  mainly  to  increased  parasympathetic,  or  vagal, 
activity  and  that  phasic  changes  during  bursts  of  rapid  eye 
movement  are  associated  with  phasic  events  in  both  the  para¬ 
sympathetic  and  sympathetic  nerves  connected  to  the  heart.  In 
other  words,  both  types  of  autonomic  activities  play  a  role  in 
physiologic  changes  during  sleep. 

There  is  a  wide  variety  of  automated  methods  of  detecting 
and  quantitating  sleep.  Larsen (38)  suggests  that  motivation 
behind  automated  sleep  stage  research  is  to  1)  provide  a  better 
understanding  of  the  physiologic  aspects  of  sleep  and  2)  to  test 
the  adequacy  of  variables  thought  to  contain  quantitative  infor¬ 
mation  about  sleep. 

Computerized  sleep  staging  involves  the  seeking  out  of 
definitive  measures  of  the  EEG  and  then  developing  some  analog. 


L2 


digital  or  hybrid  method  to  use  the  measures  for  classif ication . 
In  essence,  researchers  have  attempted  to  devise  automated 
methods  of  pattern  recognition  especially  designed  for  sleep 
patterns  of  the  EEG. 

Approaches  to  the  solution  of  the  problem  have  employed 

(23  36  42  45)  /pi \ 

techniques  in  spectral  analysis  '  '  '  ,  period  analysis1 

(28  38  60) 

and  baseline  cross-analysis v  '  '  .  These  procedures  have 

met  with  only  limited  success  except  for  the  hybrid  system 

(22) 

developed  by  Gaillard  ,  et  al. , who  reported  overall  accuracy 
(when  compared  with  visual  scoring)  up  to  92%.  An  overview  of 
the  results  of  several  studies  is  shown  in  Table  I.  Problems 
of  classifications  are  cited  as  overlap  of  measures  of  discrimi¬ 
nation^^,  the  need  of  inclusion  of  more  measures^^'^'^^  , 

(28) 

imprecise  definition  of  variables  ,  unreliability  of  visual 
(42)  .  . 

scoring  ,  and  the  need  for  additional  methods  of  analysis. 

All  of  the  above  methods  of  automated  sleep  staging  require  the 

processing  of  extensive  amounts  of  data  and  considerable  amounts 

(29) 

of  computer  time 

Our  research  program  is  dedicated  to  reducing  the  data  bulk 
and  yet  providing  accurate,  reliable,  quantitative  measures  of 
sleep.  Using  heart  rate  as  our  input  variable  has  certain 
advantages  over  previously  cited  measures.  For  example: 

1)  the  ECG  is  an  easily  measured  physiologic  parameter, 

2)  practically  speaking,  the  ECG  alleviates  the  need 
of  uncomfortable  multiple  electrodes  placed  in  the 
scalp  used  in  monitoring  EEG,  which  are  cumbersome 
and  restrictive  of  movement,  and 


TABLE  I 


REPORTED  RESULTS 
OF  AUTOMATED  CLASSIFICATION 


Stage 

Roessler^ 

T  w  42 
Lubin 

T  38 

Larsen 

Itil28 

W 

87% 

64% 

91% 

74% 

1 

45% 

70% 

64% 

73% 

REM* 

- 

56% 

66% 

- 

2 

63% 

82% 

85% 

55% 

3 

79% 

0% 

85% 

62% 

4 

84% 

91% 

85% 

59% 

OVERALL  ' 

69% 

60% 

79% 

65% 

Most  studies  did  not  attempt  to  classify  REM  separately. 


THE  PROBLEM 


The  objective  of  this  research  is  the  development  of  a 
process  by  which  rapid,  inexpensive  determination  of  level  of 
sleep  can  be  performed  accurately  using  an  easily  derived 
physiologic  parameter  such  as  the  beat-by-beat  heart  rate. 

OUR  PREVIOUS  WORK 

Our  first  study  reported  computer  sleep  stage  classification 

using  heart  rate  data  and  this  was  an  overview  of  a  two-year 

research  effort  in  the  bandwidth  reduction  of  sleep  informa- 

tion  .  Our  second  study  investigated  intrasubject 

variability  that  we  observed  during  our  first  study'  '  .  This 

study  involved  variations  of  heart  rate  during  sleep  as  a 

function  of  sleep  cycle.  A  third  study  provided  another 

integral  part  of  our  research  plan.  It  was  dedicated  to  the 

detection  of  REM  1  sleep  stage  and  eye  movement  from  beat-by- 
(58) 

beat  heart  rate  .  In  this  study  we  hypothesized  that  sudden 

cyclic  changes  in  heart  rate  observed  during  sleep  were  caused 

by  the  same  neural  process  which  caused  rapid  eye  movements. 

Our  fourth  study  involved  an  attempt  on  our  part  to  generalize 

the  work  of  the  previous  years  to  other  subjects  and  to  check 

(40  41) 

intersubject  variability  '  .  Our  classification  technique 

used  what  we  have  since  come  to  call  the  "beatquency  domain". 

(41) 

Our  fifth  study  utilized  our  previously  developed  algorithms 

on  a  wider  data  base  to  check  for  intrasubject  variability  with 
repeated  measurements  of  sleep  on  the  same  subject  over  a  7-week 


16 


period  of  time.  In  addition,  work  was  performed  using  a  super 
vised  and  an  unsupervised  pattern  classification  techniques ^ 
Our  final  project  of  this  study  was  dedicated  to  a  better  under 
standing  of  the  difference  between  awake  and  sleep  heart  rate 
data  on  a  new  data  base  that  was  collected  here  at  the 


University  of  Texas. 


METHODS 


AWAKE/SLEEP  DATA  BASE 

Analog  sleep  pattern  recordings  were  made  on  seven  normal 
subjects.  These  analog  tapes  were  recorded  in  our  sleep 
laboratory  at  the  University  of  Texas  during  August  1977  and 
contain  a  total  of  five  channels  per  subject:  one  channel  of 
ECG,  two  channels  of  EOG,  and  two  channels  of  EEG.  In  addition, 
the  14  channel  magnetic  tapes  contain  a  slow  time  code  (hrs/min) 
and  a  fast  time  code  (modulated  1  kHz  sinusoidal) .  These  data 
were  recorded  on  continuous  paper  strip  charts  which  we  used 
for  the  visual  scanning  and  scoring  by  an  experienced  sleep 
scorer.  Each  night  was  hand  scored  by  an  expert  using  the 
Manual  of  Standardized  Terminology  Technique  and  Scoring  for 
Sleep  Stages  of  Human  Subjects  by  Rechtschaf fen  and  Kales 


DIGITIZATION  OF  DATA 

Three  subjects'  analog  ECG  records  were  selected  for  digi¬ 
tizing.  For  each  subject,  comparable  amounts  of  awake  and  sleep 
data  were  digitized  as  follows: 

Record  5-36  min  of  awake,  40  min  of  sleep 

Record  6-68  min  of  awake,  66  min  of  sleep 

Record  7-32  min  of  awake,  32  min  of  sleep 

The  interval  between  successive  ECG  R-waves  (the  R-to-R 

interval)  was  measured  and  recorded  using  a  set-up  schematically 
represented  in  Figure  3.  The  ECG  analog  records  were  played 
using  a  Sangamo  14  channel  FM  tape  recorder,  amplified  to  2^ 


Figure  3 

Block  Diagram  of  Data  Gathering  System 


volts  and  processed  by  a  custom  built  interface  prior  to  entrance 
into  the  Nova  1210  digital  computer.  The  analog  to  digital 
conversion  was  performed  by  the  Nova  1210  minicomputer.  The 
output  of  the  computer  consisted  of  elapsed  time,  R-to-R 
interval,  printed  on  a  Texas  Instruments  Silent  700  output  unit 
equipped  with  digital  cassette  recording  format.  The  digital 
output  from  the  Nova  computer  consisted  of  a  string  of  R-to-R 
intervals  and  were  recorded  on  a  digital  cassette.  These  tapes 
were  subsequently  edited  by  a  interactive  digital  computer 
program  on  the  Xerox  data  system  SDS  930  computer  (Appendix  A) . 
After  these  data  were  visually  and  manually  edited  to  remove 
obvious  errors,  the  R-to-R  interval  data  were  processed  by 
program  MINDATA  into  a  format  compatible  with  our  classification 
program. 
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Introduction 


Securing  analog  heart  rate  data  free  from  noise  is 
quite  difficult  in  long  term  recording.  In  a  digitizing 
process  such  as  R-to-R  interval  recognition  this  noise 
exhibits  itself  as  sharp  spikes  and  is  often  recognized  by 
the  digitizing  program  as  R-waves.  Consequently,  it  is 
placed  in  the  digitized  data  as  small,  but  valid  data  values. 
Conversely,  infeasibly  large  data  values  may  also  appear  in 
the  digitized  data;  these  are  due  to  weak  ECG  signals  which 
are  not  detected  by  the  digitizing  program.  To  combat 
these  errors  a  method  was  needed  to  find  and  correct  invalid 
data  points  with  a  minimum  amount  of  error  added  to  the  pro¬ 
blem.  Thus,  this  program,  the  Heart  Rate  Data  Editing 
Program,  with  the  use  of  the  Xerox  930  graphics  terminal 
was  written  to  „earch  the  data  for  incorrect  values,  alert 
the  user  when  one  was  found  and  allow  him  to  make  corrections. 

The  editor  gives  the  user  several  basic  functions  in 
which  to  manipulate  the  data.  Initially,  the  user  may 
choose  what  per  cent  difference  is  acceptable  between  data 
points  and  how  many  lines  of  data  he  wishes  to  view  on  the 
screen.  With  the  aid  of  the  terminal's  function  keys, 
the  user  may  view  all  the  lines  of  data  he  has  read  in  or 
divide  the  data  into  groups  of  five  lines  only  to  avoid 
the  incessant  blinking  of  the  scope's  refresh  routine  when 
too  many  points  are  exhibited  on  the  screen  at  once.  For 
correctional  purposes,  the  user  may  divide  up  a  data  point 
which  is  too  large  or  add  points  together  which  are  too 
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small  using  the  terminal's  light  pen.  The  user  may  go  back 
and  look  at  up  to  ten  previously  processed  lines  of  data,  call 
for  more  lines  of  data  to  be  read  in,  go  to  the  next  file  or 
stop  the  program  completely. 


Tapes 

This  program  calls  for  two,  seven  track  magnetic 
tapes,  the  digitized  data  tape  for  input  and  a  tape  for 
corrected  output.  The  input  data  needs  to  be  in  card 
image  with  ten,  four  to  five  digit  numbers  per  image  (1015 
format)  recorded  at  800  BPI.  Individual  files  need  to 
be  separated  by  one  end-of-file  marker  with  the  last  file 
having  two  end-of-f iles.  The  output  data  will  also  be 
written  in  card  image,  1015  format.  A  heading,  "REC  i", 
i=l ,2, 3, . . .n,  will  preceed  each  file.  Except  for  two 
end-of-f iles  on  the  end- of  the  last  file,  no  intermediate 
end-of-f iles  will  written  on  the  tape. 


Keyboards 


Alphanumeric  Keyboard... 

The  scope's  main  keyboard  and  screen  is  basically 
that  of  any  CRT  device.  What  the  user  types  in  is  displayed 
in  the  lower  left  hand  corner  as  he  types  it  in.  If  a 
mistake  is  made  before  the  carriage  return  is  given,  the 
delete  key, "DEL",  on  the  upper  left  hand  side  will  delete 


the  entire  line  typed  in.  The  carriage  return  key,"CR", 
is  located  on  the  middle  right  hand  side  of  the  keyboard. 
Once  the  carriage  return  is  given,  what  the  user  types  in 
is  removed  from  the  screen. 


Function  Keyboard... 

To  the  right  of  the  scope's  alphanumeric  keyboard  is 
the  function  keyboard: 


<v> 

STOP 

t 

1* 

ftXT 

ft 

& 

<L 

!> 

A<X 

L 

R. 

tsc 

For  this  program  these  keys  have  been  assigned  as  follows 


Initial ization 


Upon  execution,  the  following  message  will  be  displayed 
on  the  scope. 

HEART  RATE  DATA  EDITING  PROGRAM 
PRESS  KEY  "GO"  TO  BEGIN 

The  computer  then  waits  for  a  response  from  the  user. 

Upon  pressing  the  "GO"  function  key,  the  program  will 
display 

MAXIMUM  PER  CENT  DIFFERENCE  ACCEPTABLE= 
and  then  waits  for  the  user  to  type  in  any  free  formatted 
number.  Once  it  recieves  it,  the  computer  exhibits  this 
number ( N ) 

MAXIMUM  PER  CENT  DIFFERENCE-N 
IS  THIS  CORRECT (Y  OR  N) : 

in  case  a  wrong  value  was  inadvertently  typed  in.  The 
program  will  wait  until  it  recieves  a  "Y"  or  "N".  Upon 
a  negative  response,  it  will  loop  back  and  exhibit  the 
editing  program's  title  and  wait  for  a  "GO"  function 
command.  Upon  a  yes,  it  will  store  this  per  cent  difference 
value  and  go  on  to  seek  the  number  of  lines  the  user 
wishes  th  read  in.  This  difference  value  can  not  be 
changed  without  stopping  the  program. 


Line  Request 
The  program  will  ask 


HOW  MANY  LINES  OF  DATA 

(10  LINES  MAX)  DO  YOU  WISH  TO  SEE: 


It  will  wait  for  a  number  greater  than  zero,  but  less  than 
or  equal  to  ten.  This  is  the  actual  number  of  lines  of 
data  read  off  the  data  tape.  The  data  is  not  displayed 
without  a  function  command  or  an  error  found  amongst  the 
data  values. 


Time 

An  approximation  of  the  minute  containing  the  R-to-R 
values  currently  being  displayed  is  calculated  from  the 
sum  of  all  the  data  values  read  in  and  is  written  in  the 
upper  right  hand  corner  of  the  screen  as  "i  M",  i=l,2,«..n. 


Error  Calculations 

To  find  digitizing  errors,  the  program  will  calculate 
the  per  cent  difference  between  every  two  data  points  in 
the  current  data.  If  a  difference  is  found  to  be  greater 
than  the  maximum  accepted  difference  value,  the  program  will 
quickly  try  to  determine  wjj£jch  data  point  of  the  two  is  in 
error  and  encase  it  in  a  rectangle  to  alert  the  user.  At 
times,  it  will  be  difficult  for  the  machine  to  determine 
which  is  the  incorrect  point  and,  thus,  a  point  which  is 
valid  when  considering  surrounding  points  may  be  marked 
as  an  error.  At  all  times  the  user  must  make  the  decision 
to  change  a  point,  otherwise,  marked  or  not,  the  point  will 
remain  in  the  data  as  it  is. 


If  several  errors  are  found,  the  program  will  display, 
at  most,  five  lines  of  data — two  lines  preceeding  the  line 
containing  the  first  incorrect  data  value  and  two  lines 
after  the  line  containing  the  error.  However,  if  more  than 
five  lines  were  requested,  all  the  lines  may  be  viewed 
with  errors  marked  by  giving  the  DISPLAY  command  once  the 
function  prompt  is  given  (see  below).  If  no  errors  are 
found  in  the  current  data,  the  program  will  exhibit  the 
message 

NO  ERRORS  DETECTED 

x 

and  go  on  to  request  a  function  command. 


Maximum  and  Minimum 

Every  time  all  the  current  data  is  displayed,  two 
points  xtfill  be  encased  by  diamonds.  These  are  the  minimum 
and  maximum  data  values  of  this  particular  data.  If  only 
a  fraction  of  the  requested  data  is  exhibited,  minimum  and 
maximum  points  may  not  be  marked  if  they  are  not  in  the 
few  lines  specified. 


End  of  Files 
End  of  Program 

When  an  end-of-file  is  encountered  while  reading  in 
data,  the  program  displays  the  message 
END  OF  DATA 

in  the  left  hand  corner  of  the  screen.  The  last  data  value 


is  found  and  the  program  resets  variables  accordingly  for 
the  function  commands.  The  user  is  expected  to  correct 
the  data,  if  necessary,  and  give  either  the  NEXT  or  the 
STOP  command  before  reading  in  more  data  (the  FORWARD 
command)'.  Failure  to  use  one  of  these  two  commands  results 
in  the  individual  files  read  in  to  be  tacked  together. 

If  the  NEXT  command  is  used,  and  no  other  file  exists  on 
the  input  tape,  the  second  end-of-file  is  read  in  and  the 
program  promptly  ends  with  an 
END  OF  EDITING  PROGRAM 
menage  for  the  user. 

Function  Commands 

A  prompt  for  a  function  command  is  simply 
FUNCTION: 

written  in  the  upper  left  hand  corner  of  the  screen. 

GO. .  • 

GO  is  used  only  in  the  initialization  routine  to  start 

the  program.  A  message, 

THIS  FUNCTION  CAUSES  THE  PROGRAM  TO  WAIT  FOR 
FURTHER  INSTRUCTIONS 

is  given  to  the  user. 

DISPLAY... 

DISPLAY  exhibits  all  the  lines  of  data  most  recently 
requested.  It  includes  the  computer-found  incorrect  points 
encased  by  rectangles,  the  minimum  and  maximum  values 


encased  by  diamonds  and  the  lines  containing  the  fiftieth, 

O 

hundredth,  and  hundred  and  fiftieth  elements  of  the  current 
data  marked  by  a  on  the  left  hand  side  of  the  line. 

CLEAR ... 

CLEAR  removes  everything  from  the  screen  except  the 
function  prompt,  the  number  of  the  minute  in  which  the  user 
is  currently  involved,  and  an  END  OF  DATA  message,  if 
necessary. 

THIRDS . . . 

THIRDS  allows  the  user  to  split  the  data  into  five  line 
groups  in  order  to  view  it  more  closely  without  the  continual 
blinking  caused  by  the  scope* s  refresh  routine.  The  program 
will  ask 

DO  YOU  WISH  TO  SEE  THE  FIRST  (TYPE  1), 

SECOND  (TYPE  2),  OR  THIRD  (TYPE  3)  5 -LINES 
OF  THE  DATA: 

If,  for  instance,  only  five  lines  were  originally  requested 
and  the  user  asked  to  see  the  second  or  third  five  lines 
of  the  data,  the  program  would  continue  to  display  the  above 
request  until  it  got  the  proper  response  of  one.  Otherwise, 
the  data  will  be  grouped  into  lines  1-5 »  6-10  or  lines 
11-15. 

DIVIDE. . . 

DIVIDE  splits  a  large  data  value  into  smaller  values. 

The  program  will  first  request  the  user  to 
DETECT  NUMBER  TO  BE  DIVIDED 


To  do  this,  the  user  must  hold  the  light  pen  on  the  screen, 
aim,  and  "shoot"  the  incorrect  data  value  (N  at  location  xx). 

A  "cross"  will  be  placed  on  the  value  "hit"  and  the  program 
will  do  the  following  error  checking. 

ARRAY  ELEMENT**  xx 

IS  THIS  CORRECT  (Y  OR  N) : 

The  user  must  check  the  data  value’s  location  in  the 
current  array  of  data  points.  This  is  facilitated  by 
the  " y  marks  on  the  lines  containing  elements  41-50, 

91-100,  and  141-150.  If  the  location  given  is  not  that  of 
the  desired  value,  the  program  will  ask  the  user  to 

DETECT  AGAIN 

At  times  the  computer  will  not  accept  values  "shot"  in 
the  proximity  of  the  "cross".  In  this  particular  case, 
the  computer  simply  doesn't  accept  the  "shot”  and  waits 
for  another.  Unfortunately,  it  doesn't  give  any  kind  of 
an  error  message  to  indicate  it  is  waiting.  Thus,  the  user 
needs  to  move  the  pen  down  and  over,  staying  within  the  box, 
if  necessary,  and  "shoot"  again. 

When  the  detected  value  is  the  one  desired,  the 
program  will  ask 

WHAT  DO  YOU  INTEND  TO  DIVIDE  BY  (8MAX) : 

It  waits  for  a  number  greater  than  zero,  but  less  than  eight (h) 
to  use  as  the  divisor.  However,  if  the  user  has  inadver¬ 
tantly  gotten  into  this  routine,  he  must  nonetheless  detect 
a  number,  but  can  type  in  zero  at  this  point  and  be  delivered 
back  to  the  function  prompt  without  disturbing  the  data. 


The  quotient  and,  hence,  the  new  values  (a)  will  be  calculated 
and  displayed  along  with  the  original  value  detected. 
QUOTIENT=a 

OHIGINAL  VALUE  SHOULD  BE=N 
IS  THIS  CORRECT (Y  OR  N) : 

If  this  value  is  correct,  the  array  is  expanded,  these 
values  included  and  control  returns  to  the  function  prompt. 

If  this  is  not  correct,  the  program  will  prompt  the  user 
to  type  in  the  correct  values  with 

WHAT  M  VALUES  DO  YOU  WISH  TO  USE: 
and  xfaits  for  him  to  type  M  numbers  in  one  at  a  time. 

After  all  the  numbers  are  read  in,  the  program  displays 
them. 

M  numbers 

vmTyT' 

ORIGINAL  VALUE  SHOULD  BE=  N 

ARE  THESE  THE  VALUES  YOU  WISH  TO  USE  (Y  OR  N) : 

If  these  values  are  correct,  they  are  included  in  an 
expanded  array  and  control  is  returned  to  the  function 
prompt.  If  these  are  not  correct,  the  program  loops  back 
to  the  beginning  of  this  routine  without  disturbing  the  data. 

ADD ... 

ADD  combines  numbers  which  are  too  small  to  stand  as 

individual  data  points.  The  program  will  request 

HOW  MANY  NUMBERS  DO  YOU  INTEND  TO  ADD 
TOGETHER  (10  NUMBERS  MAX) : 

and  will  wait  for  a  number  less  than  ten,  but  greater  than 
zero  to  use  as  the  number  of  addends  (M).  However,  if  the 
user  happened  into  this  routine  unintentionally,  he  can 
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leave  at  this  point  by  typing  in  zero;  this  returr£  control 
to  the  function  prompt. 

Upon  receiving  the  number  of  values  to  be  added  together, 
the  program  will  direct  the  user  to 
DETECT  M  NUMBERS* 

These  numbers  should  be  detected  in  consecutive  order 
from  left  to  right.  To  do  so,  the  light  pen  should  be 
held  next  to  or  on  the  screen,  aimed  at  the  desired  value  and 
"fired".  The  above  prompt  will  be  exhibited  on  the  screen 
until  the  correct  number  of  values  has  been  detected.  An 
error  check  will  then  be  asked  of  the  user. 

M  numbers 
NUMBERS :  'V'x'y^Z^ 

WILL  BE  ADDED  TOGETHER  TO  GIVE  A  SUM* 

IS  THIS  CORRECT  (Y  OR  N) : 

If  these  are  not  the  correct  values,  the  program  will  return 

to  the  beginning  of  this  routine.  However,  if  the  values  are 

correct  a  check  is  -also  made  to  determine  if  the  numbers 

were  detected  in  consecutive,  ascending  order  by  location. 

If  they  were  not,  the  following  message  is  displayed 

NUMBERS  NOT  DETECTED  IN  LOCATI OHALLY , 

CONSECUTIVE,  ASCENDING  ORDER 

and  control  is  again  returned  to  the  beginning  of  the  routine. 

When  the  desired  numbers  have  beu*  detected  and  added 
correctly,  the  data  array  is  condensed  to  hold  their  sum  and 
control  is  returned  to  the  function  prompt. 


*Kore  information  on  point  detection  can  be  found  in  the  divide 
command  information. 


BACK . . . 


BACK  allows  the  user  to  edit  up  to  ten  lines  of  data 
which  have  already  been  processed  for  output.  On  the  first 
read  of  a  file,  there  is  no  previous  data  to  concatenate 
with  the  present.  Thus,  the  program  will  respond  with 

NO  NUMBERS  HAVE  BEEN  PROCESSED  TO  LOOK  AT 
and  return  control  to  the  function  prompt.  Otherwise,  the 
program  will  go  on  to  ask  the  user 

HOW  MANY  LINES  BACK 

(10  LINES  MAX)  DO  YOU  WISH  TO  SEE: 

and  waits  for  a  number  greater  than  zero,  but  less  than  or 

equal  to  ten.  However,  only  fifteen  lines  of  data  can  be 

exhibited  on  the  screen  at  any  one  time.  The  maximum  number 

of  data  points  secured  from  the  preceding  data  will  be 

diminished  to  allow  the  user  to  view  all  the  current  values 

in  context  of  fewer  previous  points;  x'ive  lines  back  is 

the  maximum  number  one  can  view  with  a  full  screen  of  ten 

lines.  Also,  the  program  will  compensate  for  the  user 

asking  for  more  lines  than  are  stored.  In  this  case,  it 

will  display  all  the  lines  that  have  been  buffered  for 

output  with  the  current  data. 

When  lines  are  secured  from  the  previous  read,  they  are 
placed  in  continuity  with  the  current  lines  of  data  and  the 
entire  array  is  scanned  for  errors  as  if  all  the  data  had 
just  been  read  in.  All  functions  can  be  performed  on  this 
large  array  except  another  BACK  function.  When  a  second 
BACK  command  is  given  on  an  already  expanded  array,  all 
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the  "old"  data  points  are  replaced  in  the  memory  buffer  and 
the  original  data  array  is  restored  before  the  BACK  routine 
is  initiated.  Therefore,  the  same  data  will  be  displayed 
on  the  screen  as  before. 

FORWARD. . . 

FORWARD  stores  the  current,  corrected  data  to  allow  the 
user  to  obtain  more  data  to  edit.  It  clears  the  screen, 
prepares  the  old  buffered  data  for  output,  writes  this 
prepared  data  on  tape,  and  buffers  the  current  data.  If  the 
BACK  command  was  executed  on  the  current  data,  FORWARD  also 
restores  the  original  array  before  buffering  it. 

NEXT... 

NEXT  obtains  a  new  file  for  the  user.  -  It  first  error 
checks  to  make  sure  the  user  intended  to  execute  this  command 
by  asking 

ARE  YOU  SURE  YOU  WANT  THE  NEXT  FILE  (Y  OR  N) : 

If  this  routine  was  initiated  unintentionally,  control 
returns  to  the  function  prompt.  Otherwise,  all  remaining 
data  in  the  input  tape's  current  file  is  skipped,  data  in 

Ore. 

all  memory  buffers  as  well  as  the  current  data  xe  written 
on  the  output  tape,  all  buffers,  variables,  and  flags  are 
reinitialized  for  a  new  file  of  data  and  control  is  returned 
to  the  function  prompt.  If  this  command  was  given  in 
compliance  with  the 
END  OF  DATA 

message,  the  immediately  succeeding  file  is  the  next  one 


read  and  this  error  message  is  removed  from  the  screen. 
STOP... 

STOP  allows  the  user  to  halt  the  program.  To  protect 
the  user  from  stopping  unintentionally,  it  first  asks 
ARE  YOU  SURE  YOU  WANT  TO  STOP  (Y  OR  N) : 
and  control  is  returned  to  the  function  prompt  upon  a 
negative  response.  Otherwise,  all  previous  and  current 
data  are  written  on  the  output  tape  and  the  screen  is  cleared 
except  for  the  message 

END  OP  EDITING  PROGRAM 


Conclusion 


The  Heart  Rate  Data  Editing  Program  was  written  with 
the  intention  of  correcting  errors  in  digitized  R-to-R 
interval  heart  rr.te  data  and  was  used  to  successfully  edit 
several  files  of  data.  Programming  a  computer  to  find  almost 
all  possible  incorrect  points  and  correcting  these  points 
interactively,  increases  the  amount  of  data  wich  can  be 
edited  at  any  one  time,  decreases  the  number  of  persons 
needed  to  do  this  job,  and,  therefore,  decreases  the  amount 
of  human  error. 


'iw  n  o  c  ni\j  -  n  o  o  «-»  o  r»r»or»r»oor»  ooooooorto 
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HEART  RATE  EDTTTMG  PROGRAM 

ENABLES  USER  TO  INTERACTIVELY  EOIT  DlHITTZEiO  0 A T A 

input'. 

OIGTTTZEO  DATA  IN  1015  FORMAT  WITH  EOF  MARKINGS  BETWEEN  FILES 
OUTPUT/.  ", 

CORRECTED  DATA  IN  1015  FORMAT,  EACH  FILE  PRCEEDED  BY  SREC  XX® 

COMMON  NDATA(?00)  NPLOT  f5«T  ,  NRljF  ( 1  00  )  ,  KBUF  fl  0  0  )  ,  NTEMP  (  1  0  )  ,  MBUF  ( 1  00 

1  ) 

NDATA/:  INTEGER  ARRAY  TNTO  WHICH  DATA  IS  READ,  EACH  DATA  POINT  OCCUPYING 
ONE  WORD  IN  THE  ARRAY.  OLD  DATA  IS  WRITTEN  OVFR  UPON  EACH  READ, 

NPLOT  ’./INTEGER  ARRAY  HOLDING  THE  LOCATIONS  OF  POINTS  THE  PROGRAM  FINDS  R 
ERROR,  EACH  LOCATION  OCCUPYING  ONE  WORD  IN  THE  ARRAY,  OLD  DATA  IS 
WRITTEN  OVER  UPON  NEW  DATA  READ  IN 

COMMON/KMOD AL / IOV , MODE 
COMMON/ltm/LImI ,ltm? 

C0MMnN/F0REGN0/JCR(l?5,l  ,JSBd5PBT  , LITE  (201 
INITIALIZE  VARIABLES  AND  DATA  ARRAY 

DATA  NUM, LAST, NO/P ,0,0/NFROM, IRT1 , IPT2/ 1 HM , 0 , P/NEOF/ 1 HN/NRE AD/ 1 HN / 

DATA  LSUM,NRK,NST0P,NTTTLF,KST0P/0,O, 1HN, 1 , j  HN/NOMORE/1HN/ 

I0V=1  %  MODEsl 

NBUF(1)=F  *  KBUF(1)=B  S  MBUFC1)sO 
DO  111  .1=1,200 
11  NDATAfJlrp 

INITIALI7F  SCOPE, FORMATTED  READ  AND  WRITE  PACKAGES 
CALL  SCOPE  ( JSR  ,  1  500! 

CALL  SWRITE(IOV,0,MOOE, n 
CALL  SRFADf?) 

ASK  USER  FOR  PERCENT  DIFFERENCE  ACCEPTABLE 
CALL  INTTTALfNPOl 

•WRITE  RFC  NO’.  ON  OUTPUT  TAPJF  TO  DIFERENTI ATE  BETWEEN  FILES 
WRITE  («,2>  NTITLE 
FORMAT(*REC*, 1 X, 121 

ASK  USER  FOR  number  OF  OATA  L*f|ES  HE  WISHES  TO  SEE 
CALL  LlNFfNFROM,NLINES,NUM) 

EOFgET  SEARCHFS  FOR  EOF  AND  JUMPS  TO  EOF  ROUTINE  WHEN  IT  FINDS  EOF 
CALL  F0FSFT(6S,3) 

READ  IN  R-ROATA  POINTS 

READ(T,*M  fNDATAfKl,Ksl,NuM) 

FORMAT  C 1 0151 

NOMORFstHN  %  NREAD=IHY 
GOTO  5 

6  CALL  FOF (NUM,NDATA1 

NfOFs 1 HY 

IF  CNOMORfLeo/ HY)  GOTO  35 

C  TTME  DISPAYS  THF  APPROXIMATE  MINUTE  BEING  DISPAYED  IN  THE  UPPER  RIGHT  CORNER 
5  CALI.  TIME  (  NO  A  T  A  ,  L  SUM  ,  NUM  1 

C  FTGtIRF  FI  KIDS  DATA  POINTS  AROVF  ACCEPTED  PER  CENT  DIFFERENCE  LIMIT  AND 
C  STORES  THETR  LOCATION  IN  NPLOT 

19  CALL  FIGURE  CMDATA, NUM, NPLOT, NO, NPD, LAST, KMN, KMX, NRFAD) 

C  PROMPT  IISFR  TO  CHOSE  A  FUNCTION  KEY 

20  WRITE fj,  31) 

31  F0RMAT(1/,2X,*FIINCTTON<*) 

c 

C  USER  DEFINES  FUNCTION 

GOTO  f21 ,22,23, 2/1,25,26, 27,28,20,301  KFYM01 


c  *  •  18 

C  GO  TS  USED  ONLY  )'N  THE  INITIALIZATION  ROUTINE 
•21  WRITE  f 1 , 32) 

32  FORMAT(0/,*THIS  FUNCflON  CAUSES  THE  PROGRAM  TO  WAIT  FOR* , / , 5X , *FUR 
1THER  INSTRUCTIONS*! 

CALL  DEL  A  Y  f  3 ) 

GOTO  20 
C 

C  DISPLAY  CURRENT  DATA  ON  SCREEN 

22  LIM1  =  1  $  L I m  2  s  N  U  M 

C AH.  RFM0VEf2)  $  CALL  REMOVE  f  3 ) 

CALL  WRTTLIMfNDATA! 

CALL  PLOTTER (K mu, KMN)  *  CALL  PLOTTER  fKMX , KMX ) 

IFfwn'.EG’o)  GOTO  3R 
DO  33  J=1>0 

33  CALL  PLOTTER  f3,NPL0TU)) 

39  GOTO  ?P 

C 

C  CLEAR  screen 

23  CALL  REMOVE ( 2 )  $  CALL  REMOVE (3) 

GOTO  ?0 

C 

C  THIRDS  DIVIOFS  DATA  INTO  THIRDS  FOR  EASIER  VIEWING 
2«  CALL  REM0V£f2)  S  CALL  REMOVE  f 3! 

CALL  THIRDS  (NUM,NDATA,N0#NPL0T#KMN#KH5O 
GOTO  20 
C 

C  DIVIDE  SPLITS  A  NUMBER  WHICH  IS  TOO  BIG 

25  CALL  DIVTDE(NDATA,NUM) 

N0=0  $  GOTO  19 

C 

C  ADD  COMBINES  NUMBERS  WHICH  ARE  TOO  SMALL 

26  CALL  ADD  fNOATA,NlJM) 

NO=R  %  GOTO  19 

c 

C  RACK  UP  PROCESSED  DATA 

27  IF  fNRKltQ.cn  GOTO  UU 

IF  fNFR0M>E.1HB)  GOTO  «3 

CALL  FORRACK  fKNUM.NLOW, TPT1 , IPT2, K V AL , NUM , ND AT A , KBUF , NBUF , NTEMP , N V 
1  AL) 

«3  CALL  REM0VF(2!  $  CALL  REMOVE  (3! 

C  BACK  RECAILS  PREVIOUS  LINES  EXHIBITED 

CALL  RACKf  NDATA.'NUM,  IPT1 ,  I PT  2  ,  NBUF  ,  KBUF  ,  NTEMP  ,  KNUM ,  NLOW  ,  NV  AL  »  K  V  AL  , 
INFPOM.'NB  K> 

N0=0  $  GOTO  19 

C  MO  DATA  PPOCESSED  TO  LOOK  AT 
««  WRITE(1,«5) 

95  FORMAT (/,2X,*N0  NUMBERS  HAVE  BEEN  PROCESSED  TO  LOOK  AT*) 

CALL  DEL  AY ( 3 ) 

GOTO  ?0 

C 

C  FORWARD--GET  MOPE  DATA 
2fl  IF  (NFROmInE. 1HR)  GOTO  40 

NFROMslHM 

CALL  F0R8ACK  (KNi|M,NLOW,  TPT1  ,  IPT?#KVAL,  NUM,  NDAT  A,  KBUF,  NBUF,  NTEMP,  NV 
1AL) 

C  FORWARD  STORES  AND  OUTPUTS  CORRECTED  DATA 

90  CALL  FORWARD (MDATA, NUM, IPTt , IPT2,NRUF,KBUF, NTEMP, NSTOP,KS TOP, MBUF) 
CALL  REMOVE!?)  $  CALL  RFM0VE(3) 

IF  fMSTOPlEDl 1HY)  GOTO  39 

IF  (KSTOpIeo. JHY)  GOTO  35 

NOsB 
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•  ‘  NRK*  (NUH/Wn+NBK 
GOTO  /4 
*C 

C  GET  NEXT  FILE 
2R  CALL  OUIT  f  NSTOP,  1HF1 

XFCNSTQP'eqIihn)  goto  20 
GOTO  ?8 

3«  CAU  REMOVED 

I F ( NEOF 1 FO . 1 HN )  CALL  SKIPF (3, 1 ) 

NTITLEsNTITLE+1 
C  REINITIALIZE  variables 

M0NnRE=1HV  5  NEOFbIHN 
NFR0H=1HM  S  IPTtsO  S  IPT2SPI  S  NBK  =  0  $  NUM  =  0  $  LASTsO 
LSUMsfl  $  MBUFM)  =  0  S  KBUF(1)=0  5  MBUF(1)=0  *  NO=0  $  NSTOP=iHN 
GOTO  1 
C 

C  STOP  EDIT  PROGRAM 
30  CALL  QUTT (KSTOP, 1HS) 

IF  (KST0p1E0‘.  1  HN)  GOTO  20 
GOTO  28 

35  WRITEH.’RT 

9  FORMAT!/, *FND  OF  EDITING  PROGRAM*) 

CALL  DELAY  (5) 

STOP 

END 

C 

C 

C 

SUBROUTINE  INITIAL  (NPDl 
C 

C  INITIAL  asks  USER  FOR  ACCEPTED  PER  CENT  DIFFERENCE  LIMIT 
C 

C  input'.'. 

C  (IN  COMMON) 

C  IOVnSCREEN  OVFRLAY  NUMRFR  INDICATING  WHICH  ONE  TO  WRITE  ON 

C  MOOEsFLAG  TO  INDICATE  WHJCH  FORMATED  WRITE  MODE  TO  USE 

C  OUTPUT,’., 

C  NPD=MAXIMUM  PER  CENT  DIFFERENCE  BETWEEN  DATA  POINTS  ACCEPTABLE 

C 

CnMMON/KMODAL/IOV,MODE 

C  DISPLAY  TITLE  on  SCREEN  AND  WAIT  FOR  USER  TO  START  PROGRAM 

1  WRlTEfl,20) 

20  F0RMAT(1X#*HFART  RATE  data  EDITING  PROGR AM* , 2/ , 5X , *8  AUGUST  1R77*, 
1«/,*PRESS  KEY  *G0*  TO  BEGIN*) 

C  WAIT  FOR  USER  TO  START  PROGRAM 
N  =  KEY  (  1 ) 

WRITE (1 ,2) 

2  FORMAT  (5X,*maXIM(jm  PER  CENT  DIFFERENCE  ACCEPT  ABLE  =  *  ) 

READ  (2,3)  NPD 

3  FOPMAT(T) 

C  ERROR  CHECK 

WRITE  (1,fl)  NPD 

fl  F0RMAT(5X,*maXIMUM  PER  CENT  DIFFERENCE** , 15 ,/, 5X ,* I S  THIS  CORRECT( 

1Y  OR  N ) < * ) 

5  REAn  (2,ft)  IANS 

6  FORMATfAH 

IF  (IANS.nf'.IHm'.and'.IANS.NE^IHy)  GOTO  5 

IF  (IANS.EO’.IHN)  GOTO  1 

RETURN 

END 


ooor»oor»r>oooo  nnn  n  o  -i  lti  c  wn  w»*  o  o  o  f»  r»  o  o  r>  r»  o  <*» 
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SUBROUTINE  l  INE(nFpom,nHNES,NUM) 

L INF  FINDS  NUMBER  OF  LINES  USER  WISHES  TO  SEE 
INPUT'./. 

NFROMzFLAG  WHETHER  BACK  HAS  BEEN  CALLED  OR  NOT { B= YES , MsNO ) 
NLINESsMUMBFR  OF  LINES (  OR  NUMBER  OF  LINES  BACK)  USER  WISHES  TO  SEE 
NtJMsNUMBER  OF  ELEMENTS  TO  BE  READ  INTO  ARRAY  NDAT A 
(IN  COMMON) 

ICVsSCPEEN  OVERLAY  NUMBER  INDICATING  WHICH  ONE  TO  WRITE  ON 
MO'-EsFLAG  TO  INDICATE  WHICH  FORMATED  WRITE  MODE  TO  USE 

C'  ••'MON/KMODAL/TGV,MODE 
IF(NFROm'egIihr)  GOTO  3 
WP  :  TE  C  1  ,2) 

FORMAT(/,*HOW  MANY  LINES  OF  D  AT  A*  ,  /  ,  «>X  ,  *  (  1 0  LINES  MAX)  DO  YOU  WISH 
I  TO  SFE<*) 

GOTO  5 

OUTPUT  LINE  FOP  BACK  ROUTINE 
WRITE (1 i 4) 

FORMAT (/, *HOW  MANY  LINES  BACK*, /,5X> *( lft  LINES  MAX)  DO  YOU  WISH  TO 
1  SEES* ) 

RE  AO  (2  *  7)  NL I  NFS 
FORMATCI) 

USER  INPUT  ERROR  CHFCK 

IF((NLINEs'.le'.oIoR>LINEs’.GT/0)’,AND‘.NFRCM.EQ'.  1HB)  GOTO  3 
IF(N|.INFSlLE>>RiNLINES.GTil0)  GOTO  ) 

NUM  IS  NUMBER  OF  ELEMENST  IN  DATA  ARRAY 
NUMsNLINESMfl 
RETURN 
END 


SUBROUTINE  EOF(NUM,NOATA) 

EOF  TELLS  USER  EOF  HAS  BEEN-  FOUND  AND  SEARCHES  FOR  FIRST  ZERO  IN  ARRAY 
INPUT//. 

NUMsNUMRER  OF  ELEMENTS  USER  REQUESTED  TO  BE  IN  NDAT  A  ARRAY 
NDATAsDATA  CURRENTLY  BEING  EDITED 
(IN  COMMON) 

IOVsSCRFEn  OVERLAY  NUMBER  INDICATING  WHICH  ONE  TO  WRITE  ON 
MOOEsFLAG  TO  INDICATE  WHICH  FORMATED  WRITE  MODE  TO  USE 
OUTPUT// 

NUMsNUMRER  OF  NON-ZERO  VALID  POINTS  IN  NDAT  A 

COMMON /K MOD AL/IOV, MODE 
DIMENSION  NDATA ( 200) 

T0Vs5 

C  TELL  USFR  THERE  IS  NO  MORE  DATA 
WRITE  (J#i) 

t  F0RMAT(*EN0  of  DATA*) 

C  FIND  LAST  valid  data  value 
DO  2  KS1 , NUM 

IF  (NDATA  (K) '.EQ’.O)  NIJMsK-1  S  GOTO  3 

2  CONTINUE 

3  IOVsl 
RETURN 
END 
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SUBROUTINE  TlME(NDATA,LSUM,NUM) 

C* 

C  TIME  DETERMINES  which  minute  is  being  processed 

c 

c  INPUT, .. 

c  NDATA-ARRAY  OF  current  data  points  reing  processed 

C  LSUMsSUM  OF  ALU  DATA  POINTS  PROCESSED  SO  FAR  (NOT  CURRENT  POINTS) 

C  NUMSNUMRER  OF  ELEMENTS  IN  NDATA 

C  OUTPUT,,, 

C  LSUMsSUM  OF  ALL  DATA  POINTS  PROCESSED  SO  FAR  INCLUDING  CURRENT  DATA 

C 

COMMON/K MODAL/ I OV, MODE 
DIMENSION  NDATA(20U) 

C  ADD  TOTAL  no,  OF  SECONDS 
DO  1  J=1,NUM 

1  LSUM=NDATA(J)+LSUM 
C  DETERMINE  PROPER  MINUTE 

MIN=(L$Um/60«O0)+1 
CALL  REMOVE (R ) 
lov=a 

C  DISPLAY  minute 

WRITE(1,2)MIN 

2  FORMAT (35X, 13, IX, *M*) 

IOV=i  \ 

RETURN 

END 

C 

c 

c 

SUBROUTINE  PLOTTER(IVAR,I) 

c 

C  PLOTTER  purs  SQUARES  around  bad  DATA  POINTS  and  diamonds  AROUND  min  and  max 

c 

C  INPUT,,, 

c  Ivar=flag  indicating  plotting  diamondsc,gt;0)  or  squares(.le,0) 

C  I=LOCATlON  OF  DATA  ELEMENT  TO  BE  MARKED 

c 

DIMENSION  IXAR(5),1YAR(5) 

CALL  GETXY(I,IX,IV) 
if  (Ivar.gt:*)  goto  2 

C  PLOT  VECTORS  into  SQUARES  around  data  points 
I*AR(n  =  I*-5  $  IYAR(i)sIY«12 

I X AR ( 2 ) s I X"5  $  IYAR(2)  =  lY  +  36 

I X  AR ( 3 ) s  I  X  +  Rfl  S  IYAR(3)  =  IY  +  36 
IXaR(«)3IX+98  S  IVAR(a)=IY*l2 
IXAR(5)=Ix-5  $  IYAR(5)=IY-12 

3  CALL  SlPL0T2(3,IXAR,IYAR,5,-6) 

RETURN 

C  PLOT  VECTORS  INTO  DIAMONDS  AROUND  DATA  POINTS 
2  IXAR(5)3lX*a9  S  IYAR(5)slY-t2 

IXAR(l)  =  IXttt9  S  !YARU)sIY«12 
IXAR(2)sIX-8  $  IYAR(2)=IY+12 

IXAR(3)siX+a9  S  IYAR(3)slY+38 
tXAR(4)sI X+98  S  IYA«(a)sIY+i2 
GOTO  3 
End 

c 

c 

c 


SUP  ROUT  Inf:  FIGURF(NDATA,NI]M,NPLOT,NO,NPD,  L AS T , KMN , K M* f  NRfc  A  CO 

*c 

C  FIGURE  FINOS  POINTS  ABOVE  ACCEPTED  PER  CENT  DIFFERENCE  L I M I T  BY  COMPARING 
C  EVERY  TWO  DATA  VALUES 
C 

c  input’..'. 

C  NOATAsOATA  WHICH  NEEDS  TO  BE  EDITED 

C  NUMaNlIMPFR  OF  elements  IN  ARRAY  NOATA 

C  NPDsACCEPTABLE  per  cent  DIFFERENCE  level 

C  LAST=LAST  DATA  VALUE  OF  PREVIOUS  DATA  ARRAY 

C  (IN  COMMON! 

C  IOVsSCOPF  OVERLAY  NUMBER  INDICATING  WHICH  ONE  TO  WRITE  ON 

C  MODF=TNOICATES  WHICH  FORMATTED  WRITE  MODE  TO  USE 

C  OUTPUT , 1 i 

C  NPLOT=LOCATAIONS  OF  RAD  DATA  POINTS  JN  NDATA 

C  NOsNUMRFR  OF  ELEMENTS  IN  APRAY  NPLOT 

C  KMN  =  LOCATION  OF  MINTMijm  VALUED  DATA  POINT  IN  NDATA 

C  KMX  =  LOCATAlON  OF  MAXIMIJM  VALUED  DATA  POINT  IN  NDATA 

C  (IN  COMMON) 

C  LTMl=LOWFR  LIMIT!  LOCATION  IN  NDATA  TO  REGI  N  WRITING  DATA  VALUES  FROM 

C  LIMps UPPER  LIMIT!  LOCATION  IN  NDATA  TO  TERMINATE  WRITING  DATA  VALUES 

C 

COMMON/KMODAL/IOV.MODE 

COMMON/LIM/I.IM1  ,LIM? 

DIMENSION  NOATA(200) ,NPLOT(50) 

MIN=90000  S  MAX=D 
C  CLEAR  SCREEN 

CALL  REMOVE (2)  $  CALL  REMOVE(3) 

C  SET  LAST  TO  FIRST  DATA  ELEMENT  FOR  VERY  FIRST  ELEMENT 
IF  (LAST'.EO>)  LASTS  NDATA(l) 

C  SERACH  FOR  BAO  DATA  pniNTS 
DO  11  K  s 1 j  MU M 
TEMPlsNDATA (K) 

C  FIRST  ELEMENT  OF  array  IS  A  SPECIAL  CASE 
IF  (KlEn.t)  TEMP2sLAST  S  GOTO  10 

C  COMPARE  FIRST  DATA  POINT  IN_  CURRENT  ARRAY  TO  LAST  DATA  POINT  IN  THE 
C  PREVIOUS  ARRAY  (VERY  FIRST  ELEMENT  IS  COMPARED  WITH  ITSELF 
TEMP2  s  NDATA  (K»l ) 

C  TAKE  PER  CENT  DIFFERENCE  BETWEEN  2  POINTS 

10  PDIFsABSf ( (TEMP1-TEMP2)/TEMP2)*100,0) 

C  IF  PERCENT  DIFFERENCE  IS  GREATER  THAN  WHAT  USER  SET,  SAVE  LOCATION  OF  BAD 
C  POINT  IN  NPLOT 

IF (PDIFIgF.NPD)  CALL  SAVEBAD(K, NO, NPLOT) 

C  FIND  MINIMUM  ANO  MAXIMUM  DATA  VALUES 

11  CALL  MINMAX(K,MIN,MAX, NDATA, KMN, KMX) 

C  SAVF  LAST  DATA  POINT  IN  NDATA  ARRAY  FOR  COMPARISON 
IF  (NREAD’.EOilHY)  L  A  STsND  A  T  A  (  MUM  ) 

NREADslHN 

C  TELL  USER  if  NO  ERRORS  ARE  FOUND 

IF  (NOiFn’.O)  TOV  =  3  $  WPJTEd  ,5)  S  GOTO  6 

5  FORMAT (2/ #  *NO  ERRORS  DETECTED*) 

C  DISPLAY  DATA  WITH  BAD  PCINTS  MARKED 

CALL  WRITFRR  (NUM, NDATA, NPLOT, NO, KMN, KMX) 

6  lOVsl 
RETURN 


ouoouoouoooooooo 


SUBROUT  INF  RAVERAntKfK'O,  KlPLOT) 
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■c 

C  SAVFPAO  SAVES  LOCATIONS  of  bah  oata  points  in  array  nplot  ! 

C 

C  INPUT... 

C  KsLOCATTON  OF  POINT  IN  NO AT  A  FOUND  TO  BE  IN  ERROR 

C  NO=NUMBFR  OF  BAD  OATA  POINTS  IN  NPLOT 

C  NPLOTsLOCATTONS  OF  RAO  OATA  POINTS 

C  OUTPUT' 

C  NOsNUHBFP  OF  RAO  OATA  POINTS  IN  NPLOT 

C  NPLOTsLOCATTONS  OF  RAO  DATA  POINTS 

c 

DIMENSION  NPLOT (5n) 

JsK 

C  IF  LOCATION  OF  RAO  POINT  IS  THE  FIRST  OR  SFCONO  ELEMENT  IN  THE  ARRAY  SKIP 
C  FURTHER  TESTS  AND  STORE  THE  VALUES  LOCATION 
IF  (K.LF‘.?)  ROTO  12 
IF  (NO.LElm  ROTO  12 
IF  (NPLOT(Nn)lER.fK-n)  ROTO  13 
IF  fNPLOT (NO) LEO' (K-2))  JsK-J 

12  NOsNO+1 
NPLOT  f  NO)  sj 

13  RETURN 
END 

C  I 

c 

SUBROUTINE  WRTTFRR  ( NUM , NO A T A , NPLOT , NO , KMN , K MX ) 

FIND  AND  show  2  LINES  RFFORE  AND  2  LINES  AFTER  BAD  DATA  POINTS 
INPUT-'./. 

NDATAsAPPAY  OF  OATA  POINTS  CURRENTLY  BEING  EDITED 
NUM  =  NIJMRfR  OF  ELEMEMT5  CURRENTLY  IN  DATA  ARRAY  TO  RE  EDITED 
NPLOTsLOCATTONS  OF  RAD  DATA  POINTS  THAT  THE  PROGRAM  FOUND  IN  CURRENT  DATA 
NO= NUMBER  OF  BAD  DATA  POINTS  IN  NPLOT 
KMN  =  LOCATTON  OF  MINIMUM  DATA  POINT  VALUE  IN  NO  A  T  A 
KMX  =  LOCATION  OF  MAXIMUM  DATA  POINT  VALUF  IN  ND A T A 
(IN  COMmON)< 

IOVsSCOPE  OVERLAY  NUMBER  INDICATING  WHICH  ONE  TO  WRITE  ON 
MODEslNDICATES  WHICH  FORMATTED  WRITE  MODE  TO  USE 

LTMisLOWER  LIMITJ.  LOCATION  IN  NDATA  TO  BEGIN  WRITING  DATA  VALUES  FROM 
LlM2=UPPFR  LIMIT!  LOCATION  IN  NDATA  TO  TERMINATE  WRITING  OATA  VALUES 


COMMON/KMODAL/IOV, MODE 
COMMON/I.  IM/LIM1  ,LIM2 
DIMENSION  NDATAf2001 ,NPLOT(S0) 

C  CLEAR  SCREEN 

CALL  REMOVE  131  J  CALL  RFMOVF(2) 

C  FIND  2  LINES  BEFORE  AND  2  LINES  AFTER  RAO  DATA  POUTS 

CALI.  LlMTTS(NPLOT,NI|M,  1  ) 

C  SHOW  2  LINES  EFFORT  AND  2  LINES  AFTER  BAD  DATA  POINT 
CALI  WRIT  I.TMf  NDATA! 

C  PLOT  MIN,  MAX  POINTS  IF  THEY  ARE  WITHIN  LIMITS 

IF  (KMn'ge'lIMI ' AND’ kmnIlF.LIm?)  CALL  PLOTTER  (KMN.kmn) 

IF  (KMX, GF.LIM1 .AND. KMX. LF, LIMP!  CALL  PLOTTER  (KMX, KMX) 

C  PLOT  BAD  POINTS  WITHIN  LIMIT 
DO  3  J= 1 , NO 

IF(NPi.OTfJ)  iGE.lIMr.ANO/'PLOTfJ)  ,LE/lM2)CALL  PLOT  TE  R  ( 0 ,  NPLOT  (J  )) 
3  CONTINUE 

RETURN 
END 


,*.y  .V  .-.  .■% ,-v  .%  '  w  *.* 


oooooooooooo  o  r>  o 


S 1 1 8 R  0  U  T  T  N F  MlNMAX(K,MlNfMAX,NDATA,KMN,KHX) 
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C 

•c  MINMAX  DETERMINES  THE  MINIMUM  AND  MAXIMUM  DATA  POINTS 

C 

c  input'..'. 

C  K= ELEMENT  POINTER  IN  NOATA 

C  MTN  =  RFLATTVF  MINIMUM  DATA  VA|,UF  FROM  POINT  TO  POINT 

C  MAX=RFLATTVF  MAXIMUM  DATA  VALUE  FROM  POINT  TO  POINT 

C  NOATAsARRAY  OF  OATA  POINTS  TO  BE  FOITED 

C  OIJTPUT.'i 

c  KMNSMTMIMUM  number  AMONGST  ALL  OATA  VALUES  ENCOUNTERED 

C  KMX  =  MAXTMI|M  NIIMRER  AMONGST  ALL  OATA  VALUES  ENCOUNTERED 

C  MINeRELATIVE  MINIMUM  DATA  VALUE  1M  NDATA  FROM  POTIN  TO 

C  MAXsRFLATTVE  MAXIMUm  OATA  VALUE  FROM  POINT  TO  POINT 

C 

OIMFNSION  N0ATAI200I 

MIN=MINO(MTN,NDATA(K)I 

IF(MlN.EOlNOATAfK) IFMNsK 

MAX  =  MAXfl (max, NDATA  (K)T 

IF  (MAX'.EO.NDATA(KI  )  KMXsK 

RETURN 

ENO 


\ 

SUBROUTINE  L I M I TS ( NPLOT , NUN , K I 

LIMITS  FINDS  2  Ll^ES  BEFORE  AND  2  lines  AFTER  RAD  DATA  POINT 

input;;;  , 

NPLOTrLOCATTONS  OF  RAO  DATA  POINTS  FOUND  IN  CURRENT  DATA 
NJLIMsNIJMRER  OF  ELEMENTS  IN  NOATA 
KsLOCATTON  TO  FIND  LIMITS  BEFORE  AND  AFTER 
OUTPUT.’.; 

(IN  COMMON) 

LlMtrl.OWER  1. 1 M T T A  LOCATION  IN  NDATA  TO  BEGIN  WRITING  DATA  VALUES  FROM 
LlM?sUPPER  1. 1 M I T  A  LOCATION  IN  NDATA  TO  TERMINATE  WRITING  DATA  VALUES 

COmmON/lIm/i Im i ,LIM2 

dimension  nPL0T(5PI 

ILINESsNP|,OT(K)/lO 
LT  =  MOD(NP|_OT  (K)  10) 

if(lt;eo>)iltnes=ilines-i 

C  IIM1  IS  LOWFR  LTMIt 

LlMls(fTl.lNFS-l)*10)-R 

C  LTM2  IS  UPPFR  l T m T t 

LTM2s(ILINES+T)*10 

I F  C  l.  T  m  i  IlF>H.Imis1 

IF(LIM2;gT.NUM)LIM2sNUM 

RETURN 

END 


IN  ARRAY 
IN  ARRAY 
POINT 


I 


ooooooooooo  o  o  r>  m  r»  o  o  o  o  o  o  o  r* 
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SUBROUTINE  GETYY(T,TX,IY) 

GET  XY  DETERMINES  X,Y  COORDINATES  FROM  ARRAY  LOCATION; 

INPUT  I  LI 

IsLOCATIOM  OF  data  POINT  IN  DATA  ARRAY 
OUTPUT.II 

iXsX-COORDINiATr:  OF  DATA  POINT  ON  SCREEN 
IY=Y-C00RDINATE  OF  DATA  POINT  ON  SCREEN 

ILINFS=I/1 0 
MIJL  =  MOD(I  j  m 

IF(MUL  '.EQ.01  MUt  =  t0  S  I L  I  N  E  R  >•  I L I N  F  S  ~  1 

IFCI.FOln  TX  =  36  %  !Y=753  *  GOTO  35 

IXafMIJL*PO)-5a 

IY=753-(ILINFS*«8I 

RETURN 

END 


SUBROUTINE  WRTTLIM(N0ATA5 

SHOW  2  LINES  BEFORE  AND  2  LINES  AFTER  BAD  DATA  POINT 
input'.!  I 

NDATArARRAY  OF  CURRENT  DATA  VALUES 
(IN  COMMON) 

I0V  =  SC0PE  OVERLAY  NIJMRER  INDICATING  WHICH  ONE  TO  WRITE  ON 
MOOE=INDICATES  WHICH  FORMATTED  WRITE  MODE  TO  USE 

LIHIsLOWER  LImIT A  LOCATION  IN  NDATA  TO  BEGIN  WRITING  DATA  VALUES  FROM 
LIM?sUPPER  LIHIT4.  LOCATION  IN  NOATA  TO  TERMINATE  WRITING  DATA  VALUES 

common/ljm/limi  ,’lim? 
common/kmodal/iov.modf. 

DIMENSION  NBA T A (200) , KHAR (2) , K ARRA Y ( 1 3) 

KHAR(l)slHs  S  NCHAR=52 
CALL  REMOVE(?)  S  CALL  R E M 0 V E ( 3 ) 

C  GET  COORDINATES  OF  LOWER  LIMIT 
CALL  GETXY(LlMi#IX#IY) 

JVAl.sl.  IM  t  •  i 

JTEMP5(LIM2-LTMl)+l 

NTImESsJTEmP/10 

IF  (M0D(JTFMP,  10)  '.NE.0)  NTIMFS  =  NTIMFS  +  1  S  NSUBsMOD ( JTEMP ,  1  P )  *5 
IF  (MOD(N81J*#«)!nE.0)  NSUB  =  NSUB*(«-MOn(NSUB,«) ) 

C  GRAPH  NUMRERS  ON  SCOPE 
DO  4  K *  1 .NTIMES 

ENCODE  (52.3.KARRAY) (NDATA(JfJVAL) i J=1 , 10) 

3  FORMAT ( 1 015#  2X ) 

IF  (kIeQ'.NTIMFS)  NCHARsNStJR 

CALL  STITLF?  f2#IX,lY.18,0,KAPRAY,NCHAR#C) 

IY=TY-«A  $  IX=36  $  JVAL=JVAL+10 

«  CONTINUE 

C  DISPLAY  CHARACTER  MARKING  LINE  NO!  5 

CALL  STTTLEfP#!*, 5f»2, KHAR, 1,0) 

CALL  STTTLEfP, IP, 322, KHAR, 1,0) 

I0V  =  1 

RFTURN 

END 
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SUBROUTINE  THJ ROS ( NUM, ND AT A , NO, NPLOT , KHN , KHX) 

c  * 

C  THIRDS  EXHIBITS  only  the  third  of  the  data  the  USER  WISHES  to  see 

c 

c  INPUT, .. 

C  NUhsnUmBER  of  elements  in  ndata 
C  NDATA=ARRAY  OF  CURRENT  DATA  POINTS  BEING  EDITED 
C  NO=NUMBER  OF  BAD  POINTS  IN  NPLOT 

C  NPLOT sLOCATI ON S  OF  BAO  DATA  VALUES 

C  KMN=LOCATION  OF  DATA  POINT  WITH  THE  MINIMAL  VALUE  OF  THE  CURRENT  DATA 

C  KMX=LOCATION  OF  DATA  POINT  WITH  THE  MAXIMAL  VALUE  OF  THE  CURRENT  DATA 

C  (IN  COMMON) 

C  IOVaSCOPE  OVERLAY  NUMBER  INDICATING  WHICH  ONE  TO  WRITE  ON 
C  MODE=lNDICATES  WHICH  FORMATTED  WRITE  MODE  TO  USE 

C  LlMi^LOwER  LIMIT*  LOCATION  IN  NDATA  TO  BEGIN  WRITING  DATA  VALUFS  FROM 

C  LIM2SUPPER  LIMIT*  LOCATION  IN  NDATA  TO  TERMINATE  WRITING  DATA  VALUES 

C 

COMMON/K MODAL /I OV# MODE 
COMMON/LlM/LlMl,LlM2 
DIMENSION  NDATA(2tf0)#NPLOT(50) 

IB  WRITE(I,12) 

12  FORMAT(/,?x#*DO  YOU  WISH  TO  SEE  THE  FIRST  (TYPE  1),  * , /, 5* , *SECOND 
1  (TYPE  2),  OR  THIRD  (TYPE  3)  5  L INES* , / , 5X , *  OF  THE  DATA<* ) 

19  READ(2,17)NTHIRD 

17  FORMAT(I) 

IF  (NTHIRD.NEll  .  AND  >THI  RD  ,  NE  .  2*.  AND  ,  NTH  I RD  i  NE  ,  3  )  GOTO  18 
C  FIND  LIMITS  FOR  WHICH  THIRD  OF  DATA  BUFFER  USER  WISHES  TO  SEE 
GOTO  (13,14, 15)NTHIR0 

13  NFAR=MlN0(NUMf 50) 

LlMlsl  $  LIM2=NFAR  S  GOTO  16 
16  IF  (NUM.LE.51)  GOTO  18 

NfaR=MiN0(NUM, 100) 

L I M 1 s5 l  S  LIM2=NFAR  5  GOTO  16 

15  IF  (NUM.LE.101)  GOTO  18 
NFAR  =  MlN<3(N|jM,  150) 

LiMisim  $  lim2=nfar 

C  DISPLAY  LIMITED  DATA  ON  SCREEN 

16  CALL  WRITLIM(NDATA) 

IF  (NO.EQ.O)  GOTO  21 

C  DISPLAY  BAD  DATA  POINTS  WITHIN  LIMITS 
DO  20  Is J , NQ 

IF(NPL0T(I);GE.LIM1.AND.NPL0T(I).LE.LIM2)  CALL  PLOTTER (0,NPLOT( I)) 

20  CONTINUE 

21  IF  (Kmn.GE.I Im1.AND,KMN;lE.L2M2)  CALL  PLOTTER (KMN,KMN) 

IF  (KMX.GE.LIM1,AN0,KMX,LE,LIM2)  CALL  PLOTTER (KMX , KMX ) 

RETURN 

END 


nooooonooooo- 


SUBROUTINE  DIVIDE  (NQ A Y  A  #  MUM) 

DIVIDE  SPLITS  DATA  POINTS  WHICH  ARF  TOO  BIG 

INPUT'..: 

NDATAsAPOAY  OF  DATA  points  CURRENTLY  RFTNG  EDITED  CONTAINING  POINT  USER 
WISHES  TO  DIVIDE  INTO  S  E  V  F  R  A 1.  NUMBERS 
NUMsNUMBER  OF  ELENFNTS  IN  NDAT A 

output:: : 

MOAT  A  =  APR  AY  OF  CURRENT  DATA  POINTS  WITH  BAD  POINT  DIVIDED  INTO  SEVERAL 
NUMRFRS 

K'UMsNO,  OF  ELEMENTS  IN  NO  A  T  A  WTTH  POINT  HAVING  REFN  DIVIDED  OUT 

COMMON/FOPEGNO/JCRn  25), JSBC  J 500) , LITE f ?H) 

COMMOM/L  TM/l.iMl  f  LTN? 

DIMENSION  NO  AT  A  1200)  #L  TEMP  MR) 

WR  ITEM,  10) 

10  FORMAT  I/,*DETECT  NUMBER  TO  BE  DIVIDED*) 

C  DETECT  ARRAY  LOCATION  OF  BAD  POINT  FRM  LIGHT  PEN 

CALI  detf.ctikf.l) 

C  ASK  USER  WHAT  HE  WISHES  TO  DIVIDE  BY 
CALL  DIVNi|M(NruTA,wEL,MniV,LTEMP) 

C  DIVIDING  RY  NP6  TNDICATFS  NEED  YO  LEAVE  ROUTINE 
IF  fMOIV.FO'.O)  GOTO'  lfl 
C  EXPAND  OATA  ARRAY  TO  ENCOMPASS  NEW  POINTS 
CALL  EXPAMD(KEL,NUM,NDATA,MDIV) 

INC2=KEL 

INCtsO 

C  PUT  NEW  DATA  TNTO  DATA  ARRAY 
DO  11  Isl.’MDIV 
NOATAHNCPIsLTEMPd) 

11  INC?=INC?+t 

lfl  CALL  REMOVE  f  2  ) 

RETURN 
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m 


SUBROUTINE  niVNUH(NDATA,KPL  ,MniVf  LTFMP) 


C 

c 

c 

c 

c 

c 


DIVNUM  ASKS  USER  WHAT  HE  WISHES  TO  DIVIDE  BY 


INPUT... 

NDATASARHAY  OF  CURRENT  DATA  VALUES 
KFL=LOCATION  OF  NUMBER  TO  BE  OIVIOEO 


a 

c 

LTEMPsVALUES  TO  BE  SUBSTITUTED  FOR 

NUMBER 

'■V 

C 

OUTPUT.' 1 

C 

MOIVsDIVISOR  of  large  NUMBER 

c 

DIMENSION  NOATA(20P1 ,LTEMP( 101 

13 

NTOTAL=0 

t' 

C 

GET  WHAT  USER  WANTS  TO  DIVIDE  BY 

WRITE!!  ,  11 

1 

FORMAT!/,  1X,'*WHAT  DO  YOU  INTEND  TO 

DIVIDE 

sa 

10 

RFAO  (2,20)  MD I V 

20 

FORMAT { T 1 

IF  (MDIvIfoIbI  GOTO  6 

*  .  •  «  *  »  ^  _  r 

r: 

IF (MOIV.GT.fi, OR. MOIV.LT «P«OR,MDIV, 

Eo'.n 

C 

DETFRMINE  QUOTIENT  OF  DESIRED  NO.  DIVIDED  RY 

NOUOaNDATA (KEL1/MDIV 

GOTO  IB 

NOi  NO',  USER  DFSIRES 


12 

5 


STORE  REPLACEMENT  VALUES 
DO  3  Isl , Mpt v 
NTOTAL=NTOTAL  +NQUO 
LTEMPf  nrNOUO 

ASK  USER  IF  VALUES  ARE  ACCEPTABLE 
WRITEd  ,«)NOUO,NTOTAL 

FORMAT!/, IX, *OUOTIENT  =  *, 15,/, UNORIGINAL  VALUE  SHOULD  PE=*,I5,/,i 
10X,*IS  THIS  CORRECT  (V  OR  NlS*) 

READ!?, 51  KANS 
FORMAT (Ail 

IF(KAMS'.NE'.IHN.ANONANs'.NE.tHY)  GOTO  12 
IFfKANS'.EO.IHYl  GOTO  6 
IF  VALUES  NOT  ACCEPTABLE  ASK  FOR  VALUES 

ntotal=b 

DO  8  K=1,M0IV 
WRITE  (l,7)MOIV 

FORMAT  !/, 1X,*WHAT*,I2,*  VALUES  DO  YOU  WISH  TO  USEl*) 

RFAD(2,2)  LTEMPfK! 

FORMAT  (II 

NTOTALsWTOTAL+LTEMPfKl 
DISPLAY  ON  SCREEN  WHAT  VALUES  USF.R  ENTERED 

WRITEM.BI  MOIV,(LTEMP(J1,J  =  1,MOIV),NTOTAL 

FORMAT(/,NI5,/#*ORIOINAL  value  SHOULD  BE  =*,I5,/,*ARE  THESE  THE  V A 
1  LUES  YOU  WJSH  TO  USE  (Y  OR  N)S*) 

REAO  (2,111  LANS 
FORMAT ( A  !  ) 

C  ERROR  CHECK 

IE  (LANS.NFl  IHy'.AND'.LANS'.NE'  1HN)  GOTO  l« 

IE  (LANS.FO’.IHN)  GOTO  13 
6  RETURN 

END 


2 

6 

C 

1« 

9 


11 


aV. 


.v 
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SUBROUTTmF  DETECT  (KEL) 

PFTFCT  FINOS  LOCATION  OF 'pAA  PT  FROM  LIGHT  PEN  HIT 
OUTPUT/.  ’ 

KELSLOCATION  OF  DATA  POINT  USER  HIT  WITH  PEN 

C0MM0W/FpREGNn/.!CRM25>,  JS6Cl5P0),LITEf20) 

CALL  TRACK  (LITE, 2D) 

CALL  ONTRACKfPP) 

CALL  HIT (  T  X ,  I Y ) 

X  AND  Y  COORDINATES  PFTFRMINFP  FROM  LIGHT  PEN  HIT 
IF(  (IX-36)  >£>)  GOTO  a 
NXsl  $  ROTO  2 
NXsf  (  (1X-36) /lB)/5)+l 
NYsf (777-IY)/?a V? 

DETERMINE  LOCATION  TN  DATA  ARRAY  FROM  X*Y  COORDINATES 
KElsfNYMCM+NX 
WRITE  (1 , 3 1 KFL 

FORMAT ( / , IX,* ARRAY  ELEMENTS* , I  5 , / , 5 X , *  I S  THIS  CORRECT  (Y  OR  N)S*) 
REAPf P,6TNANS 

format(  ah 
CALL  REMOVE  ( l ) 

IF  (NANSjNE' IHvIaND'nANS.NF' JHN)  GOTO  5 
IF  (NANS.EO'.lHYl  GOTO  S 
GET  ANOTHER  PEN  IF  ELEMENT  UNACCEPTABLE 
WRITE  (1,71 

FORMAT(/,2X,*PETECT  AGAIN*) 

GOTO  1 

CALL  OFFTRACK 

RETURN 

END 


SUBROUTINE  EXPAND(KEL,NJJM,NDATA,MDIV) 

Expano  enlarges  data  array  to  hold  new  points 
input* / 

KEL=LOCATION  OF  ELEMENT  DIVIDED 
NUMsNIIMRER  PF  ELEMENTS  IN  NDAT A 

NDATAsCURRENT  data  values  along  WITH  NALUE  TO  BE  DIVIDED 
MDIVsPIVIWOR  OF  LARGE  DATA  VALUE 
OUTPUT./ 

NDATAsPATA  VALUES  WITH  ROOM  FOR  LARGE  MARKED  VALUE  TO  BE  DIVIDED 
NUMsNIIMRER  OF  ELEMENTS  IN  NEW  NPAT A 

DIMENSION  NDATA  (2PO)  ,NHOLO(20Pi) 

INClsKF.1. 

C  S A VF  ORIGINAL  DATA  ARRAY 
DO  26  Jsl(*NUH 
26  NHOLP(J)=NPATA(J) 

C  EXPAND  ARRAY  POINTER  NLIM  TO  ENABLE  DATA  ARRAY  TO  HOLD  NFW  NUMRERS 
NUMsNUM*  (MOTV-1) 

C  SHIFT  DATA  ARRAY  ELFMEMTS  TO  THE  RIGHT  TO  ENCOMPASS  NEW  NUMBERS 
DO  ID  I s  K  E  L ♦ M  0  J  V , N  U  M 
10  MOATA(nsNHPLD(I«(MDIV.l  M 

RETURN 
END 


1 

i 


.  SUBROUTINE  Ann  f  ND  AT  A  ,  NI.JM  5  . 

c 

C  ADO  COMPILES  DATA  POINTS  WHICH  ARE  TOO  SMALL 

C 

c  input!. I 

C  NOATAsCUPRENT  DATA  POINTS  WITH  THF  POINTS  USER  WISHES  TO  ADD  TOGETHER 

C  NUMeNDMRER  OE  ELEMENTS  IN  NOATA  WITH  NUMBERS  TO  BE  ADDED  TOGETHER 

C  MN  COMMON  I 

C  IOV=SCOPE  OVERLAY  NUMBER  INDICATING  WHICH  ONE  TO  WRITE  ON 

C  MOOF=TNDICATES  WHICH  FORMATTED  WRITE  MODE  TO  USE 

C  LIM1=L0WER  LIMIT!  LOCATION  IN  NO  A  T  A  TO  BEGIN  WRITING  DATA  VALUES  FROM 

C  LIM?  =  I.IPPEP  LIMIT!  LOCATION  IN  NOATA  TO  TERMINATE  WRITING  DATA  VAWUES 

C  OUTPUT.'.. 

C  NDATAsClJRRENT  DATA  VALUES  WITH  MARKED  POINTS  ADDED  TOGETHER 

C  NUMsNEW  NUMBER  OF  ELEMENTS  IN  NOATA 

C 

COMMON/FOREGND/JCBC 125} , JSBC15B0) .LITE (20) 
common/i.im/limi.lim? 

COMMON/KMODAL/TOV# MODE 
DIMENSION  NDATA  C2BB} .LSAVF  C 10} 

C  FIND  NUMBER  OF  POINTS  USER  WISHES  TO  ADD  TOGETHER 
1  CALI,  ADDNtIM  f  NADD) 

C  ZERO  NO.  INDICATES  NEED  TO  LEAVE  ROUTINE 

if  fNAon  Ieo'.  oi  goto  ?i 

C  DETERMINE  SUM  PF  THF  COMBINED  NUMBERS 

CALL  SUMfNDATA.KEL.LSAVE, KANS,NSUM#NADD) 

IF  (KAMS!eq! 1 H N )  GOTO  1 

C  ERROR  CHECK  THAT  NUMRERS  are  IN  CONSECUTIVE  ORDER 
DO  tO  KKst.NADD-l 

IF  (ISAVEfKK)  ’.LT’.LSAVE(KK  +  m  GOTO  IB 
WRITE  (t,<n  «  CALL  OEL  A  Y  C  3 }  $  GOTO  1 

9  FORMAT(/,*NUMRERS  NOT  DETECTED  IN  LOCATION  ALL Y 5V , *CONSECUT I VE 

1*  ASCENDING  ORDtR*) 

IB  CONTINUE 

C  CONDENSE  DATA  ARRAY 

CALL  C0NDENS(NADD#LSAVE,N0ATA,NL!M,KEL> 

NDATAfKFL-tNADD-n  )=NSUM 
21  CALL  REMOVE (2) 

RETURN 

END 

C 

c 

c 

SUBROUTINE  ADOnUm(NADD) 

c 

C  ADDNUM  FINDS  NUMBER  OF  POINTS  USER  WISHES  TO  ADD  TOGETHER 
C 

C  OUTPUT*.'.'. 

C  NADDsNUMBER  OF  POINTS  TO  BE  AODF.D  TOGETHER 

C 

1  wr  i  tem.  2 1 

2  FORMAT(/,*HOW  many  NUMBERS  DO  YOU  INTEND  TO  ADD* , / , 5 Y # *TOGETHER  Cl 

IB  nijmrfrs  max)<*) 

READ  f2, 31NAOO 

3  FORMA  TCI) 

if  cnadd'.eo! i !or!naod!lt,b)  goto  1 
ifcnadd!gt!ibinaoos10 

RETURN 

END 


c 

c 


cx>  -o  O'  nfc  n  o  n  -  <-n  oor>r»ooooooor' 


SUBROUTINE.  SUM(N|UTA,KEL,LSAVE,KAN8,NSU”,NAI)D) 

SUM  ADOS  BAD  DATA  POINTS’ 

INPUT  HI 

noatasappay  of  outwent  DATA  values 

NAOn  =  KiUKRER  OF  DATA  VALUES  TO  RE  ADDED  TOGETHER 
output. 1 1 

KELsLOCATTON  OF  LAST  DATA  ELEMENT  DETECTED  TO  PF  AN  ADDEND 
LSAVE=ARRAY  WITH  LOCATIONS  OF  POINTS  DETECTED  BY  USER 
KAN8  =  FLAG  INDICATING  NUMBERS  DETECTED  ARE  CORRECT  fY  OR  N) 

NSUMsSUM  OF  DETECTED  VALUES  ADDED  TOGETHER 

C0MM0N/F0RF.GND/JCB(T25T  #  JSB(15PQI  ,LITE(PP) 

DIMENSION  NDATA  f2P0>#LSAVEf  lO),MTEMPn^) 

NSUMsO 

NAODENDsNADD 

find  LOCATION  OF  EACH  RAD  DATA  PT  to  re  added 
DO  «  K=1 ,NADD 
WRI TE  C 1 » 1 D)  NAODFNO 

0  FORMAT!/, 1 Y,*DETECT*, IX. 12, 1 X , *NUM8ERS* I 
MADDENOsNAODEND-l 
CALL  DETECT  !  KEL ) 

SAVE  LOCATIONS  Or  BAD  PT 
LSAVEOO=KEL 
SAVE  VALUE  OF  RAO  PT 

MTEMPfKIsNDATA  fKEL) 

SUM  BAD  DATA  VALUES 

NSUMsNDATA  (KELUNSUM 
TFLL  USER  RESULTS  OF  SUMMING 

WRITE  f  1,  BIN  ADD,  (MTEMPOM  ,KsI,NaDD)  ,N$IJM 

FORMAT!/, *NUMBERSS*»NI5./#5X,*WILL  BE  ADDED  TOGETHER  TO  GIVE  A  SUM 
1«*#I5»/**TS  THIc  CORRECT !  Y  OR  N)S*J 
READ(2,fl)  KANS 
FORMAT! A  1 ) 

if!kans,.me'.ihn.and'.kansJ,ne‘.ihy)  GOTO  9 
RETURN 


ooooooooon  r>r>o 


subroutine  condens^add,  lsave,  ndata,  num,  keu 

CONDENS  CONDENSES  DATA  ARRAY  TO  CONFORM  WITH  COMBINING  SEVERAL  POINTS 

INPUT',.'. 

NAODsMIJ^PFR  OF  OATA  VALUES  Sl/NMfn  TOGETHER 
LSAVE  =  LOCATTON  =  LOCATinNS  OF  OATA  VALUES  SUMMED  TOGETHER 
NOATAeCIJRREMT  DATA  VALLES 
NLIMsNIIHRER  OF  ELEMENTS  IN  MDATA 
KfLsLAST  ELEMENT  DETECTED  BY  USER 
OUTPUT.^, 

NDATAsCURPENT  DATA  VALUES  WITH  POINTS  ADDED 
NIJMsNUMBER  OF  ELEMENTS  IN  ND  AT  A 


OTMFNSION  LSAVEC1PI ,NOATA(2001 
INCsKFL 

C  CONDENSE.  DATA  ARRAY 

DO  R  JrKFL-fNADD-PI  ,‘NIJM 
NDATA ( JI=NOATA (TNC+il 
9  INCsINC+l 

C  CONDENSE  NUMBER  OF  DATA  POINTS  POINTER 
NUMsNUM- f NADO-I ) 

RETURN 


cb  'i  o »  noiot  nwnw>*n  ooor>oooooor»r»ooooor>or> 
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SUBROUTINE  F0WPACK(KMUM,N1.0W,  IPTt ,  I P  T? ,  K  V  A  (. ,  N|jM .  ND  A  T  A  ,  !<  R  1.1  F  ,  NRUF  ,  NT 
1EMP.NVAL) 

FOPRACK  RESETS  RATA  ARRAVS  AFTER  HAVING  CALLED  BACK  ROUTINE 

INPUTS! 

KVAL  =  NO*.  OF  DATA  VALUES  RETURNED  to  nruf  FROM  expanded  ndata 

KMUMaNO!  OF  VALUES  TO  RE  REPLACED  INTO  KBIJF  FROM  EXPANDED  NDATA 

NLOWs LOCATION  OF  VALUE  PFFORE  1ST  NPilP  VALUE  TO  REPLACED  FROM  NDATA 

NVAL^TOTAL  no!  OF  PREVIOUS  DATA  VALUES  TO  BE  RFTIJRNFD  FROM  EXPANDED  N  D  A  T  A 1 

I P  T  1  s  N  1.1  M  MODULO  1  0  ( MUMPER  OF  ELEMENTS  IN  NT  F  MP ) 

lPT?r.NUwPFR  OF  ELEMENTS  (  IN  UNITS  OF  10)  IN  NGUF 

NUMsNUMPFR  OF  ELEMENTS  IN  EXPANDED  NDATA 

NDAT  AsARHAY  OF  CURRENT  DATA  VALUES  PLUS  SOME  PREVIOUS  VALUES 
NRUFsRUFFER  OF  CURRENT  CORRECTED  DAT*  POINTS  (UNITS  OF  10) 

NTEMPsAPRAY  OF  TRAILING  END  OF  VALUES  (MUM  MODULO  10) 

KRUFsINTFRMEDI ATE  PUFFER  (100  VALUES)  TAKEN  IMMEDIATELY  FROM  NPUF  WHEN  IT 
IS  FULL 

OUTPUT.! ; 

NDATAsDATA  values  in  their  PRE-BACK  FORM 
NUMsNUMPER  OF  ELEMENTS  IN  NDATA 

DIMENSION  KRIJF  (  1  00)  ,  NRUF  (  100)  »  NTEMP  (10).  NDAT  A  (  200  ) 

IF(KN1|m;fd>)  GOTO  3 
RESET  KDATA  ARRAY 
DO  2  I  =  1,KNIIM 

KRUF(  (100-KNIJM)  +  I)=NOATA  (I) 

CHECK  TD  SFF  IF  NBIJF  TS  EMPTY 
IFfKVAL'.LF'.O)  GOTO  5 
RESFT  NRUF  DATA  ARRAY 
DO  U  J=1 #KVAt 

NRUF  C  J  +  NLOW)=MDAT  A  (KNIJM  +  J) 

CHECK  TO  SEF  IF  NTEMP  IS  EMPTY 
if(tpti;if>)  GOTO  7 
RESET  NTEMP  ARRAY 
DO  fe  K=1 .  TPT1 

NTEMP(K)=NDATA(KNUM+KVAL+K) 

RFSET  NDATA  ARRAY 

DO  fl  Msl.MUM.NVAL 
NDATA  (M)=NDATA(NVAL<>M) 

NIJMsNUM-NVAL 

RETURN 

ENO 
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SUBROUTINE  RACK(NOATA,NUM,IPT1,IPT2,NHUF,KBUF,NTEMP,KMUM,NLOW,mVAL 
1,kval,nfrom,nbk) 

BACK  COMBINES  PREVIOUS  LINES  OF  DATA  WITH  CURRENT  DATA 

INPUT... 

ND4TA=CURRENT  data  VALUES  BEING  EXAMINED 
num=number  of  elects  IN  NDATA 

IPTlsNUH  MODULO  10  (NUMBER  OF  ELEMENTS  IN  NTEMP) 

ipt2=number  of  elements  cin  units  of  id)  in  nbuf 

NBUF-BUFFER  OF  CURRENT  CORRECTED  DATA  POINTS  (UNITS  OF  10) 
KBUFslNTERMEDIATE  BUFFER  (100  VALUES)  TAK~N  IMMEDIATELY  FROM  NBUF  WHEN  IT 
IS  FULL 

NTEMPcARRAY  of  TRAILING  END  OF  VALUES  (NUM  MODULO  10) 

NBKsNUMBER  OF  ELEMENTS  WHICH  HAVE  BEEN  PROCESSED  PREVIOUSLY 
OUTPUT... 

NDATASCURRENT  DATA  VALUES  PLUS  LINES  OF  PREVIOUS  DATA 
NUMsNUMBER  OF  ELEMENTS  IN  EXPANDED  NDATA 

knumsnumber  of  values  to  be  obtained  from  kbuf  for  expanded  ndata 

NLOW=LOCATION  OF  VALUE  BEFORE  1ST  NBUF  VALUE  TO  BE  PLACED  IN  NDATA 
NVAL=TOTAL  NUMBER  OF  PREVIOUS  DATA  VALUES  TO  BE  ADDED  TO  NDATA 
KVAL=NUMbER  OF  DATA  VALUES  TO  BE  OBTAINED  FROM  NBUF  FOR  EXPANDED  NDATA 
NFROM=FLAG  INDICATING  BACK  ROUTINE  HAS  BEEN  PROCESSED 

DIMENSION  NDATA (200) , NBUF (100),KBUF( 100) , NTEMP (10) ,  KDATA(20P) 
C0MM0N/LlM/LlMl,LlM2 
NFROMslHB 

ASK  USER  HOW  MANY  LINES  BACK  HE  WISHES  TO  SEE 
CALL  LINE(NFROM,NLINES,NVAL) 

IF  (NBK.LT.NLINES)  NVAL=NBK*10 
LIMIT=NVaL+NUM 

MAXIMUM  NO,  OF  15  LINES  OF  DATA 

IF(LIMIT,GT’.  150)NVAL=(  150-NUM).MOD(  (15P-NUM)  ,  10) 

TEMPORARlRY  STORE  DAyA  ARRAY 
DO  3  K=1,NUm 

KDATA(NVAL+K)=NDATA(K) 

IF(NVAL.GT,IPT2)  GOTO  «  “ 

KNUMs0  S  KVAL=NVAL»IPTt 
NL0WsjPT2»KVAL 
IFdPTt  ,GT,0)  GOTO  6 
GOTO  6 

KNuH=NvAL-(IPTa+IPTl) 

RETREIVE  LINES  NEEDED  FROM  KBUF ( INTERMEDI ATE  DATA  ARRAY) 

00  5  I=1,KNUM 

KDATA(I)SKBUF((100-KNUM)+I) 

NLOWa0  $  KVAL«lPT2 
IF(IPT1.LE.0)  goto  0 

Retrieve  modulo  10  numbers  bding  temporarily  stored 

Do  7  JS1.IPT1 

KDATA(KNUM+KVAL*J)SNTEMP(J) 

IF(KVAL.LE.O)  GOTO  10 

Retrieve  lines  needed  from  nbijfcstorage  array) 

DO  9  Ksl.KVAL 

9  KDATA(KNUM+K)=NBUF(NLOW*K) 

10  NUmsnVAL+NUh 

C  RETURN  COMPLETED  data  ARRAY  WITH  PREVIOUS  LINES  DESIRED  TO  NDATA 
00  11  Jsl.NUH 

11  NDATA(J)skOATA(J) 

RETURN 

END 


*  SUBPOUT  I NF  FORWARn(NDATA,NUM,  IPT1  ,  T PT2  ,  NBUF' #  KRUF  ,  NTEMP  ,  N STOP  ,  K STOP 

1  /  MBLIF ) 

C 

C  FORWARD  STORES  AND  OUTPUTS  CORRECTED  DATA 

C 

c  input;;'. 

C  noata=curpfnt  DATA  VALUES 

C  NUMsNUMRFR  OF  FI  F.MNTS  IN  N D A T A 

C  IPT1=NUH  MODULO  1«  (NUMBER  OF  ELEMENTS  IN  NTEMP) 

C  IPT?5 NUMBER  OF  ELEMENTS  (I*  UNTTS  OF  101  TN  N n U F 

c  NRUF2BUFFER  of  CURRENT  CORRECTED  DATA  POINTS  (UNTTS  OF  1«) 

C  KRUFsINTFPMEDT ATE  BUFFER  (100  VALUES)  TAKFN  IMMEDIATELY  FROM  NBUF  WHEN  ITT 

C  TS  FULL 

C  NTEMP-ARRAY  OF  TRAILING  END  OF  VALUES  (NUM  MODULO  10) 

c  NSTOP  =  F(  AR  INDICATING  TO  PRINT  OUT  KBIJF,  NBIJF,  NTEMP  (END  OF  FILE  PROCESS' 

C  KSTOPsFLAG  INDICATING  TO  PRINT  OUT  KRUF , NBUF ,  NTEMP(ENE  OF  PROFP  N ) U 

C  MBUFsMEMORV  BUFFER  FOR  OUTPUT 

c  output;;; 

C  IPTlsMUM  MODULO  10  (NIIMBFR  OF  ELEMENTS  IN  NTEMP) 

C  IPT2  =  MUMRF.R  OF  ELEMENTS  (IN  UNITS  OF  10)  IN  NBUF 

C  NRUFsRUPFFR  OF  CURRFNT  CORRECTED  DATA  POINTS  (UNITS  OF  10) 

C  KBUEsINTERMEOlATE  BUFFER  (iOt?  VALUES)  TAKEN  IMMEDIATELY  FROM  NBUF  WHEN  IT 

C  TS  FULL 

C  NTEMPsARRAY  OF  TRAILING  END  OF  VALUES  (MUM  MODULO  IB) 

C  MRUFsMEMORY  BUFFER  FOR  OUTPUT 

C 

DIMENSION  NDATA(20P1)  fNBUFdflPl)  ,NTEMP(10)  ,KBUF(1B0)  ,MBUF(1P0) 

INCP=1 

IF(IPT1  ;GT>)  GOTO  U 
NUMBRaO 

IF  (iPTP.GE'.lOO)  GOTO  8 
GOTO  9 

C  PUT  NUMBERS  MOD  IP  ;  E  M  P  0  R  A  R I L  Y  STORED  FROM  LAST  TIME  INTO  BUFFER  FOR  STORAGE 
4  CALI  0LDTFMP(NUMRR,NTF.MP,N0ATA,NBUF,IPT1,IPT2) 

IF(IPT2;lt;1PI0)  GOTO  9 

C  CORRUF  WRITFS  CORRECTED  DATA  ON  A  NEW  TAPE 

8  CALL  CORBUF  (NRUF . KBUF , MruF ) 

IPT?=B 

9  IHOLP-sNUm-miimrr 

ITEMP1=MOD(THOLD,10) 

ITEMPpsTHOLD-TTFMPt 

10  IFdTFMPp'.EO.THOLO)  GOTO  12 
DO  n  Js i , I  TEMP i 

11  NTEmP(J1sN0ATA((MHM»ITEMP1)>J) 

C  PUT  CORRECTED  ND A T A  ONTO  NBUF  ARRAY 

12  DO  13  K=INCR, TTFMP2 
NRUF(IPT2+K1=N0ATA(NUMBR+K) 
if((ipt2+k);eo;ioo)  GOTO  15 

13  CONTINUE 
IPT2=TPT?dTEMP2 

19  IPTlsITEMPl 

C  initialize  data  ARRAY 

DO  17  1=1, NUM 
17  NDAT  A (I ) sO 

C' WRITE  OUT  INTERMEDIATE  buffers  WHFN  program  STOPS  or  wen  GETTING  NEW  FILE 
IE  (NSTOP^F.lHvlAND.KSTnP.MFllHY)  GOTO  18 
WRITF(«,  1  6)  (MRUF  (K)  ,K51 , 1  0(1) 

WRITE (R, 16) (KBUFfK) ,Ksl ,100) 

IF  dPT**GTjni  WRITF(fltU)(NRUFfI),Ial,TPT2) 

IF  (IPT1.GT.0l  WRITE(R, 16) (NTEMP(T) ,lsl ,IPT1) 

16  FORMAT  (1PI5) 


ooooooor>r»ooor»  r>  o  o 


18  ♦  RETURN 

15  INCRSK+1 

*C  CORRUF  WRITFS  CORRECTED  DATA  ON  A  NEW  TAPE 

5  CALI.  CORRUF  (NBUF.KRllF.MBUF) 

6  IPT2=0 

IFf  (INCP-n’.FnlTTEMP21  GOTO  14 
IPT2=«(TNCB-l) 

GOTO  12 
END 


SUBROUTINE  OLDTFMP(NUMBRf NTEMP,NDATA,MRJF,IPT1  ,  IPT?) 

OLOTEMP  STORES  PREVIOUS  MO  1.0  NO,  INTO  STORAGE  BUFFER 

input:;: 

NTEmPsARRAY  of  TRAILING  end  of  VALUES  (NIJH  modulo  10) 

NOATAsClJPRENT  oaja  values 

IPT1SNUM  MOOULO  10  (NUMBER  OF  ELEMENTS  IN  NTFMP) 

IPTPSNUMRFR  of  ELEMENTS  (IN  UNITS  OF  IP)  IN  NRUF 

output,:: 

NIJMRRsNUMPER  OF  ELEMENTS  NFEDF  FROM  NOATA  to  make  TEN  ELEMENTS  IN  N T E M P 
NRUFsPREVInilS  DATA  VALUES  PLUS  10  NEW  VALUES  FROM  N T E M P 
1PT?=INCREMENTED  0Y  10  (10  NEW  VALUES  ADDED  TO  NRUF) 

DIMENSION  NO AT A (200)# NTFMP (10), NRUF (l 00) 

CNUMBR  =Nn;  OF  FLFmENTS  FROM  NEW  DATA  ARRAY  OBTAIN  A  LINE  OF  10  ELEMENTS 
NUMRR=10-TPT1 

C  OBTAIN  NEEDED  no:  OF  ELEMENTS  FROM  NDAT A  IN  NTFMP  ARRAY 
DO  2  Ks!,N|JMRR 

2  NTFMP(IPT1+K)=NDATA(K) 

C  STORE  NFW  LTNE  OF  NUMBERS  INTO  NBUF 
DO  3  J=1 , IP 

3  NRUF(IPT2*j)=NTFmP(,I) 

IPT2=IPT2*10 

RETURN 

END 
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♦.  * 

SUBROUTINE'  CORRUF  (  NBUF ,  KBUF  ,  HpUF  ) 

C 

C  COPRUF  WRITES  CORRECTFD  DATA  RUFFFR  ON  A  NFW  TAPP. 

C 

C  INPUT'./. 

C  NRUFsRUFFFR  of  irk  CURRENT  AND  ni.n  DATA  POINTS  TO  BF  PUT  INTO  ;<m,iF 

C  KRUFsINTFPM^niATE  DATA  BUFFER  (USED  frDR  BACK  ROUTINE)  TO  GO  INTO  MQUF 

C  NBUFsMEMORY  DATA,  RF A D Y  FOR  OUTPUT 

C  OUTPUT//. 

C  KBUFsFORMFP  MRUF  VALUES 

C  MRUFtFORMFR  KRIJF  VALUES 

C 

DIMENSION  NRUFflPiP)  .KRUFflOK)  ,MPUF(100) 

21  IFfMBUFm  ,FO>)  GOTO  17 

C  STORF  CORRECTED  DATA 

WRITE  16)  (MRI.IF(K)  ,K  =  1  ,  100) 

16  FORMAT ( 1 HT5) 

C  MBUF  IS  MEMORY  BUFFER A I T  CANNOT  BE  ACCESSED  BY  USER 

17  DO  18  Jsl.lKK 
MBUF  (J)=KBUF  (.1) 

18  KBUF  (J)sNPIjF  (J) 

22  RETURN 
END 

C 

C 

c 

SUBROUTINE  QUIT (MANS, NEIL! 

C 

C  QUIT  ASKS  USER  WHETHER  TO  STOP  PROGRAM  OR  GET  NEW  FILE 

C 

C  INPUT'//. 

C  NFILsFLAG  INDICATING  USER  WANTS  TO  STOP(=S)  OR  THE  NEXT  FILECsF) 

C  OUTPUT.'/ 

C  MANS=FLAG  INDICATING  YESf=Y)  OR  NO(=N) 

c 

IFINFTl’.EoL  1HS)  GOTO  2  _ 

C  ASK  USER  IF  HF  WANTS  ANOTHER  FILE 
WRITE  (1,1)  9  GOTO  IK 

1  FORMAT!/, SXt*ARF  YOU  SURE  YOU  WANT  THE  NEXT  F I L E ( Y  OR  N)<*) 

C  ASK  USER  IF  HP  WANTS  TO  STOP 

2  WRITE  (!,«) 

U  FORMAT  (  ',BX,*ARE  YOU  SURE  YOU  WANT  TO  STOP  (Y  OR  N)<*) 

10  READ(2,3)m/inS 

3  FORMAT ( A  1 ) 

c  error  check 

ifimans'.nf'.ihy’.and'.mans’.ne.ihn)  GOTO  10 

RETURN 

END 


u 


Input /Out put  Data  Example 


Input  Data. 


75? 

“  734 

674 

668 

7i  7 

743 

737 

726 

738 

75r> 

738 

761 

72? 

6Q6 

U 1 3 

}  114 

67? 

664 

679 

7o3 

695 

7n9 

669 

679 

675 

63? 

7?4 

74n 

7?6 

7?9 

776 

003 

746 

6?9 

66? 

676  ' 

653 

654 

645 

636 

683 

7 '3 

7o5 

f  4.  ? 

69Q 

7"7 

7  66 

7 1  5 

7?  1 

7  1  9 

673 

•  67 ' 

657 

691 

68? 

686 

8?3 

69  6 

784 

70S 

749 

7  6 1- 

755 

756 

754 

607 

7?7 

7?1 

733 

734 

77  P 

77] 

816 

8  33 

73? 

74  ** 

7?7 

7/.  1 

7 17 

74  1 

£83 

69  9 

69 1 

6r8 

684 

69  • 

69? 

6^9 

67  5 

6.77 

655 

663 

637 

68  6 

7r>2 

7?6 

69  6 

69  1 

635 

671 

1331  | 

667 

667| 

1  9?3 

7^4 

7  0C 

69  ■> 

7/.7 

748 

74  4 

7?8 

739 

7?7 

7  a.  5 

733 

7  •"  8 

693 

7  ft  6 

699 

7  ~  8 

736 

735 

7  «  0 

693 

753 

36? 

R?6 

785 

656 
1  22? 

7?6 
J7 ' 

74} 

_ 7  77 

754 
_7  ? 
71] 
679 
.~63? 
JJS^4 
7?7 
_7?6 
7«f) 
74  5 


Corrected  Output  Data. • 

'  7  ITT  "  /'/C  7 
67*  6> 


'64? 


?  67  6 

3  703 

6 — 7TT? 


703  7  03  PT 

7T*““^T“TT 
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